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FROM THE DIRECTOR 


This report should reach you approximately one month before our upcoming meeting at the 
Viscount Suites Hotel in Tucson on April 28. As usual, it begins with summary reports on continuing 
projects. These are followed by extended reports, primarily in the form of preprints of papers to be 
published in the future. 

We were informed yesterday that the DARPA program on Ultra-High Density Recording has a 
start date of today, March 19. As many of you know, it has three components — magnetic disks, 
magnetic tape and optical recording. Sixteen universities and thirteen companies are participating in the 
total program. The University of Arizona has joined with three companies and three other universities 
in the optical recording projects, which are primarily directed toward the development of "blue media." 
First-year costs at Arizona are projected at $686,346; the University will fund approximately one-third 
of this, while DARPA will fund two-thirds. This should substantially complement and expand our 
present program. 

LaserByte Corporation of Sunnyvale, California, became a Level-3 sponsor of ODSC on 
December 16, 1992. We welcome them most heartily. Our association was initiated by Frank Whitehead 
of LaserByte, who graduated from the Optical Sciences Center in 1976. He has been working on optical 
recording in the San Francisco Bay area for the past decade, first with 3M, then Verbatim and now 
LaserByte, a new company, started by Ingolf Sander, developing 3.5-inch magneto-optic drives. The 
major investor is Hyundai. We expect that Frank will attend our upcoming meeting. I hope he will be 
joined by Ingolf Sander and/or Robert Toda, the President and CEO of LaserByte. 

As many of you know, AMC’s Optical Product Division has been sold to MOST Inc., which is 
a subsidiary of Nakamichi Peripherals Inc. It now has a new name — MOST Manufacturing Inc. We 
expect that Bernard Bell of MMI will attend our meeting. I hope he will be joined by executives from 
MMI and its parent company. We are encouraging them to join ODSC, of course. 

We expect Jim Lemke of Recording Physics to be at our meeting. Several others who have been 
invited include T.C. Lee of Literal Corp., Dave Campbell of Hewlett-Packard, Tu Chen of Komag, 
Gordon Knight of Maxoptix, Archie Smith of StorageTek, and others. I expect this to be a particularly 
interesting meeting. We hope that you can attend and I look forward to hosting your visit. 
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MAGNETO-OPTIC SUPERLATTICE THIN FILMS: FABRICATION, 
STRUCTURAL AND MAGNETIC CHARACTERIZATION 

C.M. Falco, B.N. Engel, R.A. Van Leeuwen and J. Yu 

PROGRESS 

During this quarter we extended our studies to determine the electronic contribution to the 
perpendicular interface anisotropy in Co-based multilayers. Using in situ Kerr effect measurements, we 
investigated the influences of different transition metals (TM = Ag, Au, Cu, and Pd) on the magnetic 
properties of single-crystal Co films grown on Pd (1 1 1) and Au (1 1 1) surfaces. We reported last quarter 
the discovery of a large peak in the perpendicular anisotropy when approximately one monolayer of Cu 
or Ag is deposited on the Co surface. 

We now have added a computer-controlled stepper-motor drive to our MBE sample transfer 
mechanism. The motor allows us to move the sample at a constant velocity from behind a shutter during 
deposition. The film, therefore, is deposited as a wedge with a linear variation of thickness across the 
substrate. In this way, we can study a continuous range of coverage on a single sample. The stepper 
motor also provides the necessary control for precisely positioning the sample in the laser beam for Kerr 
effect measurements at the different coverages. 

Also during this quarter we have continued our work on MnSb alloy films. We have made 
sputtered Mn x Sb,. x (x * .50) alloy films using DC triode system with a 3N5 Mn target and 5N Sb target 
both bonded to Cu backing plates.’ The typical sputtering power for Sb was 60 Watts (15 Watt/in 2 ) and 
120 Watts (30 Watts/in 2 ) for Mn. Typical system pressure before sputtering was * 1 x IB 7 Torr. 
Normal argon sputtering pressure was 7 mTorr, although samples were made in pressures ranging from 
3 mTorr to 20 mTorr. Substrate height varied from 3.63 in. to 7.06 in. Samples were deposited on 
substrates either at room temperature or -250 °C. Substrates were normally at ground potential but in 
some instances also were biased to +55 VDC and to -55 VDC during deposition. 

Films of pure Sb and Mn approximately 1000 A thick were sputtered on carbon and 90° 
orientated A1 2 0 3 . Precise measurements of the deposition rates were determined by Rutherford 
backscattering spectrometry (+ 3%) and the film structure was analyzed by x-ray diffraction. Standard 
Bragg-Brentano 6-26 x-ray diffraction techniques showed the Sb film to be rhombahedral and 
poly crystal line with the (001) orientation dominant, although many other orientations were present. 
Bragg-Brentano spectra of the Mn film consisted of only one small peak possibly attributable to Mn. 
Subsequent measurements were performed using a Seemann-Bohlin x-ray diffractometer with a glancing 
angle of incidence of 6° so that the effective film thickness is increased to approximately ten times the 


*3N5 means 3 nines 5. It is the way buyers and sellers of materials specify chemical purity. Ed. 
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actual thickness. This technique showed that the film was polycrystalline a-Mn phase with predominant 
(330, 41 1) orientation. 

The Mn^Sb,., films were deposited by rotating the substrates over the sputtering guns and 
sequentially depositing thin layers of each of the two materials. The sputtering power was set so that the 
ratio of the thicknesses, t^/t^, was 2.47. This results in 50 at. % of each material( x = .50). With this 
setup, arbitrary layer thicknesses can be deposited by simply changing the sample rotation speed. Typical 
layer thicknesses for Mn was 2.6 A and for Sb was 6.4 A. Total films thicknesses were approximately 
1000 A. We have used different substrates to try to influence the growth of the Mn x Sb!. x . These include: 
glass, Si(l 1 1), Si(001), GaAs(l 10), A1 2 0 3 , mica and quartz. To determine the effect of different buffer 
layers on the M^Sb,., films, we deposited polycrystalline films of Mn, Sb and Cu on these same 
substrates. 

Bragg-Brentano x-ray diffraction performed on the MnSb films showed that the films were 
polycrystalline with two orientations normal to the film plane: (101) and (110). Even though the peak 
intensities were tow, their widths were narrow. Calculations of the crystallite size from the peak widths 
using the Scherrer equation yielded values between 300 A and 400 A. We did not observe any peaks due 
to Mn 2 Sb, Mn or Sb. While the single crystal substrates generally resulted in higher intensity peaks, 
possibly due to a smoother starting surface, there was no significant difference in the orientation of the 
films due to different substrates. Films grown on buffer layers of polycrystalline Sb, Mn, and Cu also 
were not significantly different from each other or from those without buffers. To date, we have 
systematically varied sputtering rates, sputtering pressure, substrate height, substrate bias, buffer layers, 
amorphous- or single-crystal substrates, and layer thicknesses. We have not yet found a technique that 
results in c-axis oriented films. 

Ex situ annealing of a film at 300 °C for 5 minutes or 400 °C for 10 minutes in an Ar atmosphere 
did not change the structure; however the film appearance degraded significantly for the higher 
temperature. A comparison of x-ray spectra showed little change except that the x-ray intensity increased 
for the MnSb(lOl) and decreased for the MnSb(103) peak after annealing. 

VSM measurements indicate in-plane magnetization with values of M, = 2.0 - 2.2 ^ B /Mn atom. 
This is lower than the bulk value of M„ which is 3.3 n B !Mn atom. Samples deposited on Si had values 
of 2.2 /r B /Mn atom, while those samples on glass were 2.0 fi B l Mn atom. In-plane coercive fields ranged 
from 100 Oe to 200 Oe. Measurements of polar Kerr rotation yielded typical values of 0.2° regardless 
of substrate type, while the polar ellipticity, e k , was much smaller at —0.02°. 

We have recently begun deposition of MnSb by MBE. The work is in preliminary stages and 
results will be reported next quarter. 
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PLANS 

During the next quarter we will extend our studies of Co interface anisotropy and the Mn/Sb 
system. We are interested in the magneto-optic and magnetic properties of c-axis oriented MnSb alloy 
films. Therefore, we will use MBE to produce highly textured films for study. From these high-quality 
films, we will be able to estimate the fundamental performance limits of the oriented MnSb alloy and, 
therefore, guide our sputtering studies of more practical media. 
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CHARACTERIZATION FACILITY FOR MAGNETO-OPTIC MEDIA 
AND SYSTEMS 

M. Mansuripur, H. Fu, S. Gadetsky, S. Sugaya, T.H. Wu, J. Zambuto, R. Gerber, 
T. Goodman and J.K. Erwin 

OBJECTIVES 

Objectives of this research are listed below. 

Measure the hysteresis loop, Kerr rotation angle, anisotropy energy profile. Hall voltage and 
magnetoresistance of thin-film magneto-optic media using our loop-tracer. This instrument has 
a maximum field capability of 21 kOe, and can measure the temperature-dependence of the 
media’s magnetic properties in the range 77 K to 475 K. 

Measure the wavelength-dependence of the Kerr rotation angle, 6 k , and ellipticity, e t , for thin-film 
media using our magneto-optic Kerr spectrometer (MOKS). The spectrometer operates in the 
wavelength range 350 nm to 1050 nm. 

Measure the dielectric tensor of thin-film and multilayer samples using our variable-angle 
magneto-optic ellipsometer (VAMOE). This device measures (at the HeNe wavelength of 
633 nm) the complex reflection coefficients of the samples at angles of incidence ranging from 
26° to 88°. A computer program then estimates the dielectric tensor from the measured 
coefficients of reflection. 

Measure the hysteresis loop, coercivity, remanent magnetization, saturation magnetization, and 
anisotropy energy constant for thin film magnetic media using vibrating sample magnetometry. 
The temperature range of our VSM is 77 K to 1000 K, and it has a maximum magnetic field 
capability of 12 kOe. 

Observe small magnetic domains and investigate their interaction with defects using magnetic 
force microscopy. 

Perform static read/write/erase experiments on thin-film magneto-optic media using our static test 
station. The test station is a polarized-light microscope, modified to include a laser diode, 
electromagnet, temperature-controlled stage, and a TV camera, all under computer control. 
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Integrate the existing models of magnetization, magneto-optic effects, coercivity, and anisotropy 
in an interactive and user-friendly environment. Analyze the characterization data obtained in 
the various experiments, using this modeling package. 

Measure focusing- and tracking-error signals on a static testbed; determine the "feedthrough" for 
various focusing schemes, investigate the effects of polarization and birefringence, and compare 
the results with diffraction-based calculations. The testbed has two HeNe lasers (X, = 544 nm, 
Xj = 633 nm ), a series of high-NA objectives (0.5 to 0.8), piezoelectric actuators (with better 
than 0. 1 nm positioning accuracy) for driving the sample, CCD camera emulating multi-element 
detectors, several interchangeable focus-error modules, and in-built shearing and Twyman-Green 
interferometers for ascertaining beam quality and/or calibrating the focus position. 

Measure the birefringence of optical disk substrates using two variable angle ellipsometers. The 
first system operates at X = 780 nm and has the capability of both reflection and transmission 
measurements. The beam in this system has an angle of incidence 6 in the range of (0, 70°), and 
an azimuth of incidence <j> which can vary continuously from 0 to 360°. The second system 
operates at X = 633 nm, is restricted to reflection measurements only, and has an angle of 
incidence in the range of (20°, 70°). The azimuth of incidence in this system is restricted to two 
values: 0 and 90°. 

PROGRESS 

Substrate Birefringence Measurements 

A significant fraction of our efforts during the past quarter was devoted to measurements of 
substrate birefringence. We now have three systems and four people devoted to this project. A complete 
report on this work appears as Appendix A and soon will be submitted to Applied Optics. The following 
is a brief description of the three systems with emphasis on their main features. 

VAE633. This variable-angle ellipsometer operates at the HeNe wavelength of 633 nm, can make 
reflection measurements only, and is limited in its setting of the azimuth of incidence (0 = 0 and 0 
= 90° are the only possibilities). The system originally was built in our laboratory for dielectric tensor 
measurements for magneto-optic media, and recently was modified for birefringence measurements. 

VAE780. This variable-angle ellipsometer operates at X = 780 nm, has both the reflection and 
transmission capability, and allows arbitrary settings of the azimuth of incidence (0 = 0- 180°). The 
system has been designed specifically for measuring birefringence on 3.5 in. and 5.25 in. optical disk 
substrates. 
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Commercial ellipsometer. This variable-angle and variable-wavelength ellipsometer was donated 
to us by Eastman Kodak Company (we thank Dr. Joseph Wrobel of Kodak Research Laboratories for 
suggesting and facilitating the donation). We have made a minor modification to the system and have 
begun measuring substrate birefringence at various wavelengths (X = 400 nm - 800 nm). This work is 
in the preliminary testing stage and we hope to present its results at the IAB meeting in April. 

Ring-Tone Lens for Focus-Error-Detection 

Work on the ring-toric lens continued during this quarter. We have fabricated several more 
lenses and tested them in our static testbed. A report on the progress of this work is included as 
Appendix B. The upgraded testbed can measure at two wavelengths (red and green HeNe) as well as two 
polarization directions (parallel and perpendicular to grooves). We have done extensive comparisons 
between the ring lens scheme of focus error detection and the conventional method of using an astigmatic 
lens. 

Magnetic and Magneto-Optic Characterization of Magneto-Optic Media 

Hong Fu has completed a series of studies on thermomagnetic write/erase cycle using the 
Connection Machine simulations. He will present an invited paper on this subject at the Intermag meeting 
in Stockholm, April 1993. He is working on a long paper for Computers in Physics based on these 
studies and expects to have the paper ready in time for the next quarterly report. 

Sergei Gadetsky, our new research scientist, arrived in December and is engaged in a study of 
Co/Pd superlattice films produced by Dr. T. Suzuki of IBM Almaden Research Center. His findings are 
summarized in Appendix C. 

PLANS 

During our studies of substrate birefringence we have noticed certain curious phenomena related 
to the relative orientation of incident beam and grooves on the substrate. We plan to look into these 
phenomena carefully and will try to understand their origins. 

The ring-lens studies will continue and we will try to make lenses with smaller ring diameters. 
We also will investigate the dependence of track-crossing signal and feedthrough on the relative 

orientation of polarization vector and the groove. 

With the graduation of Bruce Bernacki and his departure for Oak Ridge National Laboratory, we 
will miss one of our highly productive members. New students will be learning to use DIFFRACT and 
will continue Bemacki’s modeling of the optical path in magneto-optic systems. 
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MeV ION-BEAM ANALYSIS OF OPTICAL DATA STORAGE FILMS ~ 

v J.A. Leavitt, L.C. McIntyre Jr. and Z. Lin 

OBJECTIVES 

Our objectives are threefold: 1) to accurately characterize optical data storage films by MeV 
w ion-beam analysis (IBA) for ODSC collaborators; 2) to develop new and/or improved analysis techniques; 

and 3) to expand the capabilities of the IBA facility itself. 

Using 'H + , *He + and 15 N +t ion beams in the 1.5 MeV to 10 MeV energy range from a 5.5 MV 
Van de Graaff accelerator, we determine film thickness (in atoms/cm 2 ), stoichiometry, impurity 
w concentration profiles and crystalline structure by Rutherford backscattering (RBS), high-energy 

backscattering, channeling, nuclear reaction analysis (NRA) and proton induced X-ray emission (PIXE). 

Most of these techniques are discussed in detail in the ODSC Annual Report (February 17, 1987), p. 74. 

The PIXE technique is briefly discussed in the ODSC Annual Report (March 15, 1991), p. 23. 
v< 

PROGRESS 

Film Characterization 

From March 1, 1992, to March 1, 1993, we provided a total of 214 equivalent backscattering 
v analyses of thin films provided by ODSC collaborators. Table 1 indicates the distribution. 



Period 

Falco 

3/1/92 - 6/1/92 

10 

6/1/92 - 9/1/92 

27 

9/1/92 - 12/1/92 

31 

12/1/92 - 3/1/93 

19 

3/1/92 - 3/1/93 

87 


Table 1. Number of analyses of thin films for ODSC. 


Zelinski Armstron 


22 


50 


8 



Totals 


34 


98 


45 


37 


214 


For Falco (with Engel, Van Leeuwen, Wiedmann and Yu), we provided channeling analyses of 
Pd/Co multilayers, backscattering analyses of PdCu, MnSb, NbCuMn, PdCoAu, and CoAglr films and 
PIXE analyses of several films containing adjacent mass elements (such as Pd and Ag). 

For Mansuripur (with Wu and Shieh), we provided backscattering analyses of 3-layer (and 4- 
layer) SiN/TbFeCo/SiN/(AI) films to determine stoichiometry and "thickness." 

For Zelinski (with Vogt and Weisenbach), we performed channeling studies of Si0 2 /Ti0 2 (Zn) 
waveguide films to provide stoichiometry, "thickness" and impurities. 

For Armstrong (with England), we provided two backscattering analyses of SnS 2 films. 
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LEAVITT, MCINTYRE AND LIN 


Technique Development 

We have continued to develop techniques for increasing sensitivity for quantifying light elements 
in/on heavy matrices. In this connection: 1) We have published measured 4 He- 10 B, 4 He-"B 

backscattering cross sections for incident “He energies 1 MeV to 3 MeV (see publications). 2) We have 
published details of a technique for quantifying 10 B and U B in films using ( a,p ) nuclear reactions for 
incident 4 He energies near 3 MeV (see publications). 3) Work involving use of ( a,p ) reactions to 
quantify ,4 N has been submitted for publication (see preprint. Appendix A, ODSC Quarterly Report, 
December 15, 1992), 4) Work involving use of 4 He and 'H non-Rutherford backscattering to quantify 
9 Be in thin films will be reported at the 1 1th IBA Conference in Hungary, July 1993. 5) We are 
continuing work on using (a ,p) reactions to quantify A1 and Si in films. 6) We have recently developed 
a technique to accurately determine the stoichiometry of SiC films on single crystal Si substrates. The 
technique uses channeling with incident 3776 keV 4 He ions to reduce the Si substrate signal and enhance 
the C signal. The 4 He-C cross section is about seven times Rutherford for this 4 He energy. The 
hydrogen content of the films is determined by the elastic recoil detection (ERD) technique. 7) The PIXE 
technique is now fully developed. It has been used to analyze more than 2500 samples during the past 
year. Most of these samples are of biological interest; we have used the PIXE technique on 10 films for 
ODSC to separate adjacent-mass film components. 

Facility Development 

We are continuing to work on installing the nuclear microprobe. The goal of this project is to 
produce analysis beams (*H or 4 He) of about 1 jtm diameter. This microbeam will be rastered over the 
sample to determine the lateral distribution of trace elements in features more than a few microns in size. 
The beamline is nearly complete; we have actually run a beam through the entire length of the apparatus. 
We have not yet finished installing two high-vacuum pumps. An extended period of testing will follow. 
We hope to be in operation in late spring 1993. This project is a collaborative effort involving Q. 
Fernando (University of Arizona, Department of Chemistry) and is partially supported by funds from an 
NIH Superfund grant (see the ODSC Quarterly Report, September 15, 1997, p. 27). 

PLANS 

During the next year, we plan: 1) to continue characterizing optical data storage films for ODSC 
collaborators; 2) to continue our program for developing techniques for quantifying light elements in/on 
heavy matrices; and 3) to complete construction and testing of the nuclear microprobe. 
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APPLICATION OF HOLOGRAPHIC OPTICAL ELEMENTS TO 
MAGNETO-OPTIC READ/WRITE HEADS 


R.K. Kostuk, E. Campbell and T. Kim 

OBJECTIVES 

Objectives of this research are to determine the theoretical and practical performance limits of 
holographic optical elements (HOEs) formed in different recording materials, and to evaluate the 
application of these components to magneto-optic read/write heads. 

PROGRESS 

Hologram Characterization, Modeling, and Design 

We completed measurements of the refractive index profile of dichromated gelatin emulsions for 
recording volume holograms. Characterizing the refractive index profile is important for realizing 
accurate designs using our diffraction models. Measurement techniques based on Brewster angle and 
Fresnel reflection coefficients were used to determine the emulsion’s surface index. An interferometric 
method with oil index matching was used to determine the bulk refractive index of the emulsion. The 
experiments showed that the surface index was higher than the bulk index, which contradicts the change 
expected due to exposure, dehydration, and the resulting airation of the emulsion. We suspect that the 
emulsion surface collapses during development, forming a layer with higher refractive index. We found 
some references to experimental data that support this hypothesis, but we would like to confirm this 
model on our own with additional experiments and SEM analysis. A short proceedings paper 
summarizing these findings is reproduced in Appendix D. 

Gene Campbell is examining the application of volume gratings for polarization and phase 
components. His model, based on rigorous coupled wave theory, allows the index and permittivity 
modulation to be varied as a function of emulsion depth. This characterization is necessary for an 
accurate description of the phase and polarization properties of diffracted fields. Figure 1 shows the 
calculated phase difference between s-polarized and p-polarized diffracted fields as a function of 
permittivity modulation for a zero-order beam passing through the grating. As indicated, a permittivity 
modulation of about 0.17 (An ~ 0.0623) provides a X/4 phase delay, and a modulation of 0.27 
(An - 0.10) produces X/2 retardation. Both of these modulation values can be achieved in DCG, and 
it may be possible to reduce the required index modulation value by forming thicker emulsions. Figure 
2 shows the efficiency of a polarization beam splitter (PBS) as a function of grating strength with an 
interbeam angle of 60° within the emulsion. The grating strength is defined as the product of Aed/X. 
To achieve a grating strength of 1.2 to make the PBS with X = 0.50 fim and d = 8.3 pm, an index 
modulation of 0.027 is needed, which can be realized in DCG. 
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We also examined the effect of saturation on the efficiency properties of volume gratings during 
this period. The model used for a saturated refractive index modulation is 

n(7) - JV [1 -exp(-£(7)/E.)] , 


where n( r) is the saturated index modulation, N 0 is a proportionality constant, E( r) is the grating 
exposure, and E e is a saturation constant. (This model also was used by workers at Solymar’s group at 
Oxford.) Figure 3 shows the efficiency with a saturation constant of 1000 mJ/cm 2 , and N 0 of 0.5. Some 
of the important consequences of saturation is that the TE efficiency does not go to zero, and the 
non-oscillatory efficiency function for the TM polarized field. The more commonly quoted results of 
Kogelnik have TE and TM polarized fields that have a sin 2 dependence on increasing modulation and 
thickness. This difference will greatly influence the design of many types of polarization beam splitters. 


TE - TM Phase Difference 



Permittivity Modulation 
Figure 1. Zero-order holographic waveplate. 
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TE and TM Efficiencies versus Grating Strength 



Grating Strength 

Figure 2. Sixty-degree polarization beam splitter. 


60 Degreee Polarization Beam Splitter 



Exposure (mJ/cm2) 

Figure 3. Effect of permittivity modulation saturation. 
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PLANS 

Our future efforts will be directed toward: 1) designing and fabricating phase control components 
in DCG and surface relief grating elements; 2) helping a new student become familiar with the surface 
relief grating model developed by Charles Haggans (now at 3M); 3) investigating hybrid diffractive 
optical components that combine surface relief and volume type gratings; and 4) extending and verifying 
new modeling efforts that include diffraction of finite beams, converging beams, and cascaded grating 
systems. We also would like to expand our modeling capability to incorporate finite element analysis 
techniques, which are important for chirped grating analysis. 

We are currently preparing a final report summarizing our grating model work for IBM. The 
future of this program, however, is unclear since it is no longer supported by ODSC or IBM. We hope 
to secure funding for this project from other sources that are primarily directed toward optical 
interconnect applications and that also require the use of diffraction models. 
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DYNAMIC TESTBED LABORATORY AND MICRO-OPTICS 

T. Milster, K. Erwin, F. Froehlich, J. Kann, W. Li, W. Schlichting, R. Shetty, E. Walker and 
M. Wang 

OBJECTIVES 

Objectives of this research are as follows: 

Understand the three classes of superresolution and their behavior in optical data storage systems; 

Investigate new and improved components and techniques in the optical system for data detection 
and servo control; 

Apply micro-optic components to the optical system to reduce size and weight; 

Investigate techniques, such as near-field optical probes, for recording data densities beyond that 
possible with superresolution; and 

Understand and find solutions for problems associated with dynamic testing, especially those that 
arise when evaluating blue-sensitive media. 

PROGRESS 

Our group has been very busy in several areas. In this progress report, we highlight some of the 
accomplishments of this quarter. In particular, we include an extended discussion of the superresolution 
work with the AMC testbed. We also discuss plans for our near-field optical microscope (NSOM) 
instrument and some theory that describes the NSOM transfer function. The MCD work is proceeding 
nicely, our static experiments are complete and a paper is ready for submission. 

Theory 

Our theoretical work this quarter consisted of continuing to develop a simple model for the 
NSOM work. We are completing a paper that will be submitted to Optics Letters on this subject. 

Modeling 

In the last several ODSC reports, we have described the effects of superresolution by using a 
shading band in the collection aperture. Our results indicated that both signal and noise are reduced with 
the shading band, producing no significant change in the carrier-to-noise ratio (CNR) except at high 
frequencies. However, the transfer function of the system is changed so that we may be able to take 
advantage of equalizing the frequency response of the system. A summary of our results has been 
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accepted for presentation at the ISOM ’93 meeting this summer. Appendix E is a copy of the extended 
abstract. 

We are now studying a modification of the simple shading band described in previous ODSC 
reports. Our configuration is shown in Fig. 1. We shape the shading band into an "X" that is designed 
to pass the + first diffracted orders at a specific frequency. In Fig. 1, the shading band is designed for 
around 4 MHz. The portion of the beam blocked by the shading band only contributes noise. The net 
effect is an increase in the CNR at this frequency. At lower frequencies, some degradation in signal is 
observed, but the CNR is improved due to the reduced noise. At higher frequencies, the signal is not 
affected, but the noise is reduced, again improving the CNR. An additional benefit of the X-band 
shading is that the transfer function of the system can be somewhat equalized. Our scalar diffraction 
model was used to determine the effect of X-band shading on the signal. The results are plotted in Fig. 
2, which indicates a slight improvement in the transfer function with the shading band. In the Dynamic 
Testing section below we discuss our experimental results. 



Figure 1. X-band shading. The exit pupil of the optical system used to read out data is shown with the 
± first diffracted orders of the signal. The darkened portion of the picture corresponds to the shading band. 
This shading band is designed to pass most of the signal power at 4 MHz and block excess noise. At lower 
frequencies, the shading band blocks some of the signal, but the noise is reduced more, resulting in a net 
CNR increase. At high frequencies, the signal is not blocked, and the shading band only reduces noise. For 
our dynamic experiments, we tested the shading band at different offsets in the pupil. 
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Frequency (MHz) 

Figure 2. Diffraction model results with and without the X-band shading. A slight improvement in 

the transfer function is observed when X-band shading is used. 

We presented a paper at OPTCON ’92 titled "Figures of Merit for Laser Spot Quality," which 
describes the effects of aberrations on a simple gaussian laser beam that is passed through a diffracting 
aperture. An interesting observation of Ed Walker’s modeling work is that caution should be used when 
implementing the full-width-at-half-maximum (FWHM) criterion for beam quality because it is very 
insensitive to aberration. A preprint of the paper appears as Appendix F of this report. 

Static Experiments 

We are continuing our research program in near-field scanning optical microscopy (NSOM). 
Fred Froehlich is constructing a sophisticated NSOM instrument to study the application of near-field 
technology to high-density optical storage. NSOM has been used to write and read MO domains as small 
as 60 nm (see E. Betzig et al., "Near-field magneto-optics and high density data storage," Appl. Phys. 
Lett. 61(2), 1992). Our interest is to analytically model and experimentally characterize the detection 
of MO domain edges with a near-field tapered fiber-optic probe. WYKO instruments has provided an 
atomic force microscope (AFM) that we are adapting for NSOM use. Froehlich has submitted the 
abstract (Appendix G) that will be presented at the ISOM conference this summer. 
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Figure 3. Contour plots of ly for the quadrilayer structure, with the off-diagonal elements of the 

tensor restored: a) magnetization up; b) magnetization down. 

Mark Wang is studying some aspects of superresolution by examining the amplitude and phase 
of the signal component reflected from MO media. Masud Mansuripur has predicted (Appl. Opt. 30(22), 
3154, 1991) that high signal values exist in outer regions of the four quadrants in the pupil, as shown in 
Fig. 3. Besides trying to experimentally measure this effect, we will try to see the effects of moving the 
readout spot over a domain edge. The experimental aspects of this problem have proven difficult in the 
past, but we hope to improve the sensitivity by using a modification of a readout scheme that takes the 
difference signal optically (Milster and Yao, US Patent No. 4,823,220). 

Dynamic Testing 

We have configured the AMC testbed for optical equalization experiments, as reported in the 
December 15, 1992, ODSC report. Weiquan Li and Wang have been performing the tests. Our latest 
tests are with X-band shading. Figure 4 shows the results for a 3M disk at 3 MHz signal frequency. 
The graph is plotted versus offset of the shading band (see Fig. 1 for an explanation of offset). At 1.0 
mm offset from the center, the CNR increases by about 3 dBm. Improvements in CNR also are observed 
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at other frequencies. The reason that the maximum CNR does not occur in the center of the aperture is 
that the beam is not centered properly in the testbed. We will report more details on our experiments 
at the IAB meeting in April. 

Wolfgang Schl itching has received an MO disk from Kodak that is optimized for MCD detection. 
He is evaluating the medium on the dynamic testbed. Since much more light gets back to the detectors 
with MCD, there are some modifications that need to be made to the electronics to handle the additional 
power. For example, the focus servo amplifiers saturate during writing and erasing, so their gains must 
be reduced. Without this adjustment, we can get about 41 dB CNR at low frequencies. More details will 
be presented at the IAB meeting. 

In other dynamic testing news, Wang will be presenting his work on the differential wax-wane 
focus servo at ISOM ’93. (A preprint of his abstract appears as Appendix H). We also are beginning 
work on the blue-laser head that will be used with the next set of superresolution experiments. These 
experiments use the air-bearing teststand. 
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Figure 4. Carrier and noise properties versus X-band shading offset at 3 MHz signal frequency. 
Notice that the CNR is improved by 3 dB at an offset of 1.0 mm. The offset is necessary due to a 
slightly decentered beam in the exit pupil of the drive. 
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K.A. Bates, J.K. Erwin, N. Ramanujam, L. Li and JJ. Burke 

OBJECTIVES 

Our immediate objective is to understand the limitations of guided-wave and grating coupler 
devices in their application to optical data storage. Our long-range goal is to develop and validate design 
codes for integrated optic devices. 

PROGRESS 

As a member of the short wavelength laser source team (SWAT) our principal research activity 
has been in the development of numerical models for the design of a blue wavelength integrated optical 
source for data storage applications. Nanda Ramanujam’s research efforts during the last quarter have 
been focused on designing a polarization-selective Bragg reflector for an extended cavity semiconductor 
laser using a grating-assisted slab waveguide. This laser acts as a pump for quasi-phasematched second 
harmonic generation (QPM-SHG) in an indiffused-KTP segmented waveguide. The Bragg reflecting 
grating is able to effectively discriminate between the TE and TM modal fields and reflects only the TM 
field, thereby permitting only TM laser operation. Ramanujam is currently exploring designs (such as 
chirped grating periods) to further increase the reflection bandwidth. During the next quarter, he plans 
to develop models and design programs for grating-assisted devices on channel waveguides. 

Keith Bates continues to work on the variable groove depth grating coupler project. In the past 
year or so he has made grating couplers that produced near-Gaussian output beams. The fabrication 
method was based on an idea first suggested by Lifeng Li, i.e., moving a slit with variable dwelling time 
in front of the resist grating under ion beam bombardment. However, certain imperfections of the stepper 
motor that was used to drive the slit resulted in some minor intensity fluctuations in the output coupled 
beam. At the beginning of this year, Kevin Erwin suggested a new method (the invention disclosure is 
in preparation) to produce the desired variable depth grating. Initial testing of this new idea has shown 
very promising results. Bates is currently pursuing it. 

Despite the advances made during the past 25 years in the electromagnetic theory and numerical 
modeling of diffraction gratings, certain grating problems remain difficult to solve even on a super- 
computer. For example, the problem of plane wave diffraction by a highly conducting metallic grating 
of arbitrary shape in TM (electric vector perpendicular to the grating grooves) polarization cannot be 
analyzed by any existing method when the ratio of groove depth to period is greater than 2. Many 
methods develop numerical instability long before the ratio reaches 2. 

Recently, Li has constructed a computer program that is based on the classical modal method, 
its extension to the conical incidence case, and the R-matrix propagation method, which he borrowed 
from the chemical physics literature. The program approximates an arbitrary grating profile by dividing 
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it into many layers of rectangular gratings. Because of the use of the /?-matrix propagation algorithm the 
program is immune from the usual numerical instabilities associated with the presence of evanescent and 
anti-evanescent waves in the conventional propagation algorithm. In addition, because of the excellent 
convergence properties of the classical modal method, Li’s new method also exhibits good convergence 
characteristics when applied to deep metallic gratings in TM polarization. For example, the program can 
predict diffraction efficiencies of a gold grating of sinusoidal profile and groove depth to period ratio of 
100 at a near-infrared wavelength (we realize that such groove depth is not technologically possible.) 
The program can treat diffraction gratings of arbitrary surface profile, arbitrary refractive index gradients 
in the direction perpendicular to the grating plane, and arbitrary angle and polarization of the incident 
plane wave. The program could easily be converted to treat waveguide gratings and, in particular, the 
segmented planar waveguides of non-rectangular segments and non-uniform index distribution. 

This work is being summarized for submission to a journal. The next step will be to extend the model 
to include anisotropy of the grating material. 
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COMPACT DISC ERROR MEASUREMENTS 


D. Howe, K. Harriman and B. Tehranchi 

OBJECTIVES 

The objectives of this project are as follows: 

Provide hardware and software that will perform simple, real-time, high resolution (single-byte) 
measurement of the error burst and good data gap statistics seen by a photoCD player read 
channel when recorded CD write-once discs of variable quality (i.e., condition) are being read 
(9/93); 

Extend the above system to enable measurement of the hard decision (i.e., 1-bit error flags) and 
soft decision (i.e., 2-bit error flags) decoding information that is produced/used by the Cross- 
Interleaved-Reed-Solomon-Code (CIRC) block decoder employed in the photoCD player read 
channel (1/94); 

Construct a model that uses data obtained via the systems described above to produce meaningful 
estimates of output error rates (due to both uncorrected ECC words and misdecoded ECC words) 
when a CD disc having specific (measured) error statistics is read (completion date to be 
determined); and 

Check the hypothesis that current adaptive CIRC block decoders are optimized for pressed 
(DAD/ROM) CD discs. If warranted, do a conceptual design of an adaptive CIRC decoder that 
is optimized for write-once CD discs (completion date to be determined). 

BACKGROUND 

Equipment that determines whether a given CD disc conformans with the Phillips "Red Book" 
error specifications is commercially available. To make this determination, however, this equipment does 
not monitor data read from the disc, i.e., it does not directly measure disc errors. Instead, it monitors 
pin outputs on its CIRC block decoder chip (i.e., the VLSI circuit that performs the decoding of the Cl 
and C2 ECC words) that indicate when either a one-byte, or two-byte error is corrected, or when an error 
of length greater than two-bytes is detected, by either the Cl or the C2 ECC decoder. The equipment 
averages the results of this pin monitoring activity over one-second intervals before reporting them. 
Although this is adequate for a "yes/no" determination of the disc’s conformance to the Red Book error 
specifications, it is not adequate for high resolution (single byte) mapping of a disc’s errors. To see that 
this is true, consider that at the nominal CD playback scan velocity of 1.3 m/s, one second represents, 
on average, five full spiral turns of the 120 mm diameter disc, i.e., a surface area of 2.08 mm 2 (since 
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the nominal track-to-track spacing on a CD disc is 1,6 ). This should be compared to the disc real 

estate occupied by a single byte; one byte occupies a length of 5.1 /*m along the disc’s spiral track, which 
corresponds to a surface area of 8. 1 6x 1 O' 6 mm 2 on the disc. The coarseness of error event measurements 
delivered by commercial CD disc evaluation equipment, together with the fact that such equipment does 
not directly "look at" the data being read from the disc, diminishes the usefulness of this commercial 
equipment for certain applications. Activities such as disc manufacturing process development/control 
and disc aging/degradation studies require the observation of small changes in the type and number of 
errors residing on a disc. 

The system identified in the first of the four project goals listed above, which will be completed 
this summer, will provide high resolution (single-byte) measurement of the error events on specially 
recorded CD discs (i.e., discs that have a repetitive 24-byte data pattern written to them using the CD- 
Audio format) that are played on a Phillips Model CCD-461 CD-ROM player. The fact that a particular 
data pattern must be present on the disc is not a limitation since this system is meant for use with (i.e., 
studies of) write-once CD discs which can have any specified data written on them for experimental 
purposes. (Note; we have determined that our hardware could be adapted to high resolution measurement 
of errors on CD discs containing arbitrary data, e.g., pressed CD-ROM discs.) 

PROGRESS 

To understand the following description of the CD error measurement system that we are 
developing, one must know how the various CD data structures, such as Cl and C2 codewords and EFM 
frames, are organized on (i.e., recorded on) a CD disc. This background information is contained in 
Appendix I ("Determining the Locations of the Various CIRC Recording Format Information Blocks 
[User data blocks, C2 & Cl word and EFM frames] on a Recorded Compact Disc”). It is recommended 
that those interested in the details of our CD error measurement system read that Appendix to facilitate 
their understanding of what follows. 

The scheme we are using to perform high resolution measurement of the errors on CD discs is 
enabled by the architecture of the Phillips CD3A CIRC block decoder, which is a VLSI circuit that 
performs data detection, channel data demodulation and Cl & C2 ECC decoding. It is the latest 
generation CIRC block decoder; it features adaptive data detection, intelligent channel clock 
synchronization, and soft decision ECC decoding. This chip’s Cl ECC decoder will determine how to 
handle a given received Cl codeword based on the calculated snydrome for that word and whether any 
1-bit erasure flags are attached to any of the bytes of the word (by the eight-to-fourteen channel data 
demodulator); its C2 decoder chooses how to treat a given received C2 word based on the calculated 
syndrome for the word and the distribution of 2-bit soft-error flags that accompany the word (a 2-bit flag 
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is attached to every C2 word byte when the byte is output by the Cl decoder). The CD3A chip is used 
in the Phillips CCD-461 CD-ROM player as well in the Kodak PhotoCD player. 

A high level diagram of the CD3A chip is given as Fig. 1. The aspect of this chip’s architecture 
that we exploit is the fact that data are clocked from the chip itself to an external RAM, and back, as 
illustrated in Fig. 1. This is done to implement the de-interleaving that is required to obtain Cl 
codewords from the EFM frames that are sequentially read from the disc, and to obtain C2 codewords 
from the data blocks that are output by the Cl decoder. By monitoring the data traffic marked as write 
1" in Fig. 1 we intercept the serial data stream recovered from the disc after eight-to-fourteen channel 
data demodulation has occurred, but before any ECC decoding. After appropriate synchronization is 
obtained, this data stream is organized into individual EFM frames. 

Phillips CD3A Decoder 


Decoder Chip 



Figure 1. 

If a given 24-byte block of user data is repetitively recorded on the disc being read, the lower 
32 bytes of each of the recovered EFM frames will be identical, assuming that there are no errors in the 
block (cf.. Appendix I). Therefore, by comparing each recovered EFM frame, byte by byte, with the 
exact 32-byte pattern that would be present in each EFM frame in the absence of errors, any error that 
is contained in an EFM frame can be located. Fig. 2 illustrates this; the block labeled "error detector" 
in that figure represents the 32-byte comparator. This logic element outputs 32 bits in response to each 
EFM frame comparison operation; a "zero" is output if two compared bytes are identical, a "one is 
output if they are not. The hardware will be developed in two phases (versions). Phase I will produce 
one output byte for each EFM frame in which an error is detected and the output byte will specify the 
number of erroneous bytes in the contaminated EFM frame. The Phase I hardware will provide 
intermediate resolution error measurements, i.e., error counts will be averaged over a spatial disc track 
segment that corresponds to a single EFM frame, which nominally occupies a distance along the track 
of about 0.176 mm (this distance includes EFM frame overhead). An advantage of the Phase I system 







HOWE, HARRIMAN AND TEHRANCHI 


EFM 


Error 


Fram# 


Detector 


CltD 


32 

bytes 


o 

• 32 

® I 


\ 



Phase II 
4 

Bytes 


\ 


counter 


one 

byte 


Phase I 


Figure 2. 

is that its output data can be processed in real time and we plan to provide polar graph visualizations of 
the measured error statistics. Phase II hardware, which will produce 4 output bytes for each EFM frame 
comparison, will provide high resolution (single-byte) error burst and good data gap length measurements. 
Output data from this system will be stored for post processing. 

As an extension to this project , we plan to extend the Phase II hardware’s functionality to the 
collection of the 1-bit erasure flags that are generated by the eight-to-fourteen demodulator in response 
to channel data run-length violations (cf., the second project goal listed at the beginning of this report). 
These flags are passed across the CD3A CIRC block decoder RAM interface along with the EFM frame 
data (cf.. Fig. 3). We also will study the possibility of collecting the 2-bit soft decision flags that are 
generated by the Cl decoder and passed to the C2 decoder for its use in selecting an error handling 
strategy. These 2-bit flags appear on the CD3A RAM interface at the point labeled "write 2" in Fig. 1. 

During the past quarter we have determined the sequencing used for writing and reading both data 
and error flags to/from the CD3A RAM; we have devised a method of obtaining EFM frame block 
synch; we have worked out a scheme for relating measured error events to a geometrical location on the 
disc; and we have completed a high (logical block) level design of the CD error measurement hardware. 
We are currently designing/constructing the functional system subblocks and writing the system control 
software. 
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Phase II Hardware 
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Abstract 


The birefringence of bare and coated substrates for magneto-optical recording is experimentally 
investigated using ellipsometry at the wavelengths of 632.8 nm and 780 nm. We measure the 
rotation and ellipticity of the polarization state of the reflected or transmitted light for different 
incident angles and different orientations of the incident linear polarization. The measured data 
is then fitted by a computer program which solves the Maxwell equations for the plane-wave 
propagation in a multilayer structure and minimizes the error between the measured and 
calculated data by adjusting the unknown parameters of the multilayer. This approach enables 
us to determine orientations of the three principal axes in the substrate and the corresponding 
refractive indices. A special feature of our ellipsometers is that a glass hemisphere is placed in 
contact with the substrate, which eliminates the refraction of the incident light and enables a 
maximum propagation angle of 70° (with respect to the normal) in the substrate. This increases 
the sensitivity of the measurement. Certain anomalies have been observed, which we believe are 
associated with the presence of grooves on these substrates. 


t Visiting scientist from Functional Devices Research Laboratories, NEC Corporation, Japan. 
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I. INTRODUCTION 


Birefringence is a phenomenon observed in optically anisotropic media in which the refractive 
indices depend on the polarization direction. The polarization-dependent phase-shifts due to 
birefringence of the plastic substrates of magneto-optical (MO) disks can easily degrade the 
performance of MO readout systems [1-7], As the MO data storage systems progress towards 
high resolution and ultrafast speeds, the birefringence has become an increasingly important 
problem. 

Birefringence in plastic substrates is due to preferential molecular orientation and internal 
stresses. The refractive index is higher when the electric field of the light beam is parallel to the 
molecular chain direction. For injection- molded polycarbonate substrate the fluid shear during 
the injection molding process acts to align the polycarbonate chains in the plane of the disk. 
Therefore, the refractive indices n, and n 2 for the in-plane electric field components E, and Ej 
of light are much larger than n 3 for the perpendicular component E 3 . The characterization of 
substrate birefringence has been the subject of several papers [ 6 , 8 - 11 ], which confirmed the 
above picture and found that n, - n 3 is on the order of IB 4 , while | n, - n 2 1 is on the order of 
1C 5 . However, the previous work is limited in the following aspects: First, the formula [8,11] 
by which the experimental data were processed applies only to transparent single-layer media 
under the assumption that the principal axes of the index ellipsoid are fixed in the radial, track 
and perpendicular directions. The phase-shifts due to the multilayers and due to the tilt of the 
index ellipsoid could not be taken into account. Second, due to refraction at the interface 
between air and substrate, the propagation angle of the light in the substrate is usually small. For 
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example, 70° incident angle in air corresponds to 36.5° propagation angle in the polycarbonate 
substrate (n = 1.58). So the beam possesses only a small Ej component and the experimental 
sensitivity to n 3 is limited. 

The present paper describes measurements of the substrate birefringence using one ellipsometer 
with wavelength X = 632.8 nm and another with X = 780 nm. We remove the first restriction 
of the previous work by using a computer program which solves the Maxwell equations for a 
plane- wave propagating in a multilayer structure including the substrate [12]. This approach 
enables us to determine the three principal axes of the substrate and the corresponding refractive 
indices. To overcome the second limitation, we have placed a hemispherical glass lens atop the 
substrate to eliminate refraction. A maximum propagation angle of 70° in the substrate has thus 
been achieved. We measured the rotation angle and ellipticity of the polarization state of the 
reflected and/or transmitted light as functions of the incident angle for linearly polarized incident 
beam. The investigated substrate materials include glass, polycarbonate plastic (bare and coated), 
and amorphous polyolephin plastic (bare). Measurements were performed on regions with or 
without grooves. For the former case, the grooves were either filled with index-matching fluid 
or left empty. Different angles between the plane of incidence and the direction of grooves were 
used in the measurements. This paper is organized as follows. The two ellipsometer systems and 
measurement procedures are described in Sec. II, along with a description of the data analysis 
procedure using the multilayer program. The measurement results are described in Sec. III. Sec. 
IV contains some concluding remarks. 


n. ELLIPSOMETER SYSTEMS AND THE MEASUREMENT METHOD 
H.l The VAE633 System 

Figure 1 shows the experimental setup. The light source is a HeNe laser with X = 632.8 nm. 
As the beam first passes through the polarizer and the quarter- wave plate (QWP), whose fast 
axis is at 45° to the transmission axis of the polarizer, it becomes circularly polarized. The 
second polarizer selects the linear polarization direction of the incident beam. A hemispherical 
glass lens is placed in contact (using index-matching fluid) with the surface of the sample. This 
hemisphere is actually slightly less than one half of a sphere, and is mounted in such a way that 
the geometric center of the spherical surface is on the back side of the substrate. The light enters 
normally into the hemisphere, goes through the substrate, and is reflected from the back side of 
the substrate. After reflection, the beam passes through another QWP with an adjustable fast-axis 
orientation. The beam then goes through a Wollaston prism which can be rotated about the 
beam. Finally, the two beams of mutually perpendicular polarization are detected by the 
photodetectors. 

The hemispherical lens eliminates refraction at the air-substrate interface, thus producing a larger 
propagation angle within the substrate. For example, in polycarbonate substrate with n = 1.58, 
without the glass hemisphere one can reach a maximum angle of 39° at grazing incidence (i.e. 
0«c = 90°), while with the hemisphere in place (n^. = 1.52), the propagation angle inside 
the substrates is 64° at 0^ = 70°. 


H.2 The VAE780 System 


The setup is schematically shown in Fig. 2 (a). The operation principle is similar to the VAE633 
system. One of the main differences is that all the moving parts are mounted in one arm, which 
can be rotated in both polar and azimuthal directions. The alignment of the system is thus easier 
to maintain. The freedom of two-dimensional rotation is particularly valuable to measure the 
effect of grooves, which are found to be sensitive to both polar incident angle 0^ and the 
azimuthal incident angle relative to the radial direction. The linearly polarized light from 
the laser diode (X = 780 nm) passes through the beam splitter and enters normally through the 
hemisphere before it reaches the sample. The direction of incident polarization is controlled by 
rotating the laser in its socket. The hemisphere is used in different ways, depending on whether 
the sample is transparent or reflective. For transparent samples, two hemispheres are used, see 
Fig.2 (b). The one on the top is transparent and is used to eliminate refraction. The one at the 
bottom is coated with aluminum and reflects the beam normally so that the light goes back along 
the incident path. In this configuration the light passes twice through the substrate. For coated 
samples, we only need the top hemisphere, see Fig.2 (c). Half of this hemisphere is coated with 
aluminum. The beam enters from the uncoated side, gets reflected from the bottom of the 
sample, bounces off the coated side of the hemisphere, and retraces its path. Thus the beam 
passes four times through the substrate. In both transmission and reflection modes, the use of 
hemisphere(s) does not introduce phase shift or relative amplitude change between p and s 
components upon entering and exiting the substrate. Therefore, the measured values of rotation 
and ellipticity are solely determined by the properties of the sample. As the light beam returns, 
it is routed by the beam splitter into the detector channel, where its polarization state is 
determined by the use of the QWP, a Wollaston prism and two detectors. 


H.3 The Measurement Principle 

The measured quantities are the rotation angle B and ellipticity angle e of the reflected beam. 
The definitions of B and e and the corresponding measurement procedures are described below. 
The p- and s-polarization directions for incident and reflected beams in VAE630 and VAE780 
systems are defined in Fig. 1 and Fig. 2 (a), respectively. The polarization angle of the 
incident beam is measured from the p-direction, Fig. 3 (a), and can be read from the setting of 
the polarizer in VAE633 and the laser diode in VAE780. In general, due to phase retardation 
between the p- and s-components, the reflected beam is either left- or right-handed elliptically 
polarized (LEP or REP). Figure 3 (b) shows that the ellipse can be characterized by the 
orientation of its major axis, B u and the ellipticity angle e. By definition | tan £ | is the ratio 
of the minor to major axis, and £ < 0 for LEP and e > 0 for REP. Thus, -45° ^ ^ 45°. 

In particular, £ = -45° (or +45°) corresponds to LCP (or RCP), and £ = 0 corresponds to 
linearly polarized light. When the fast-axis of the QWP is parallel to the major axis of the 
ellipse, the beam passing through it becomes linearly polarized. Its polarization angle 0 W can be 
smaller or greater than B^, depending on whether the beam before the QWP is LEP or REP. 
These two cases are shown in Figs. 3 (c) and (d). B^ and 0 W can be measured simultaneously 
by adjusting the orientations of QWP and the Wollaston prism, until one of the detectors 
receives no light. Once extinction is achieved, the orientation of the fast-axis of the QWP gives 
while the Wollaston orientation gives 0 W . The rotation and ellipticity angles of the reflected 
beam are then given by: 

® - ®<ni» ~ ^ pol > c = ~ ®«i» * W 

In practice, the error for the measured value of £ is less than 1°, as estimated from repeated 
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measurements. The error for 8 depends on the value of e. For | e | < 40°, error of less than 
1° can be achieved. For | e | >40°, the ellipse is close to a circle and the orientation of the 
major axis is hard to measure. The repeatability of 8 in such cases is within several degrees. 


0.4 Data Fitting with the Multilayer Program 

The first part of the multilayer program [12] solves the Maxwell equations for a plane wave 
obliquely incident upon a multilayer structure consisting of a given number of layers. Each layer 
is characterized by its dielectric tensor and thickness. With adequate restrictions on the tensor 
elements, different types of media (isotropic, birefringent, magneto-optic, etc.) can be described. 
Figure 4(a) shows the quadrilayer structure of a typical MO disk. This structure is used in data 
fitting for coated substrates. For bare substrates, we consider the substrate alone without 
additional layers. The hemisphere material (n^j^ = 1.52) and the index-matching fluid (n^ = 
1.52) are treated as the medium of incidence. Since the optical paths in VAE633 and that in the 
transmission and reflection modes of VAE780 are different, three versions of the multilayer 
program were used in the data analysis. The second part of the program searches the best values 
of unknown parameters by minimizing the error, defined below, between the measured and the 
calculated data: 


error 
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In this equation, the superscript "exp" stands for experimental data, and "cal" for calculated 
data, N is the total number of the measurements on a given sample. The unknown parameters 
are the refractive indices, the orientation angles of the principal axes of the ellipsoid of 


birefringence and the structure parameters of the multilayer. In our Convex computer system, 
it takes only a few minutes to search for 10 unknown parameters in fitting 120 data points. 

In Fig. 4 (a) X and Y denote the in-plane radial and track directions of an optical disk, and Z the 
normal direction. In addition to the incident angle 0*, we use to specify the azimuth of the 
incident k-vector. For all the samples, we found that one principal axis (3) is always along the 
perpendicular direction, but the in-plane axes (1 and 2) may be rotated in the X-Y plane, see 
Fig.4 (b). To be specific, we let 1 be the in-plane axis with the larger refractive index, and use 
^m-pime to denote the in-plane rotation of axis 1 from the radial direction. 

There are two aspects of the actual systems that are not taken into account in the multilayer 
program. One is that the beam is not a plane wave, but a focused beam with 0.5 mm diameter 
that converges to a 50#*m spot at the bottom of the substrate. For thick samples, the two 
reflected beams from the front and back surfaces are spatially separate and do not interfere with 
each other, while in the program the two beams overlap. However, since the reflection at the 
interface between the hemisphere and the substrate is weak, the interference only causes small 
fluctuation in the calculated curves. Therefore, the plane wave approximation is acceptable. 
Another aspect is the structure of the grooves on the surface of the substrate, which are usually 
0.5 /im wide, 70 nm deep, and have a periodicity of 1.6 /tm. The effect of the grooves on 
polarized light has not been fully investigated. We found that, for bare substrates, the data 
measured at = 0° can be matched very well, while for coated substrates, good match can 
be obtained for = 90°. For other values of there is usually mismatch between the 
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measured data and theory. In these occasions we are probably witnessing anomalous behavior 
caused by the grooves. 

m. RESULTS 

The measured samples include a glass microscope slide, two 3.5" bare polycarbonate substrates 
(PC3.5B01 1 and PC3.5B02), two 5.25" polycarbonate bare substrates (PC5.25B01 and 
PC5.25B02), one 5.25" amorphous polyolephin substrate (APO5.25B01), two 3.5" MO disks 
(PC3.5C01 and PC3.5C02), three 5.25" MO disks (PC5.25C01, PC5.25C02, PC5.25C03). All 
coated disk substrates are polycarbonate plastic. The main results of these measurements are 
summarized in Table 1. 

m.l The Glass Sample 

To check the VAE633 system, we measured a glass microscope slide with thickness d = 1.2 
mm. Glass is believed to be isotropic with a refractive index close to that of the hemisphere and 
index-matching fluid. So the optical path is simple (e.g. no multiple reflections) and the problem 
should be well described by the well-known formulas for reflectivity at a single interface 
[13,14]. 


Fig. 5 shows 0 and e as functions of the incident angle in the range 10° ^ ^ 70° for ^ 

= 45°. The measured data in (a) and (b) are same. The calculated curves in (a) are obtained by 


‘In this notation PC stands for polycarbonate (APO stands for amorphous polyolephin), the number 
3.5 for the disk diameter in inches, B for bare substrate (C for coated substrate), and 01 is the serial 
number. 



assuming that the glass is isotropic. The match gives n gU „ = 1.54. We see that the rotation 
angles can be matched very well, but the non-zero ellipticity for 9 mc <41° (critical angle of 
total internal reflection) cannot be matched. The average fitting error is 1.9°. To obtain a better 
match, we allowed the glass to have a small amount of birefringence. The curves in Fig. 5(b) 
show the best match with an average error of 0.8°. The fitted refractive indices are n, = 1.54, 
(n, - n 2 ) = 6x 10 -6 , (n, - n 3 ) = -lx 10^ and ^-pUne = -9.2°. In this particular case where the 
sample has no grooves, X and Y stand for in-plane directions along the two rectangular sides 
of the microscope slide. 

The basic features of the plots of 6 and e versus 0^ can be understood using the reflectivities 
for Ep and E, for dielectric interface [13,14]. Below the critical angle for the total internal 
reflection (0^ < 41°), the reflectivities are real and non-equal, i.e. the values of | Ep j and 
| E, | will change after reflection but the phase difference should remain zero. Therefore, the 
reflected beam is still linearly polarized (e = 0) but in a different direction (0 * 0). The 
measured non-zero e indicates that there exists a small amount of birefringence. At the critical 
angle the behavior of 0 and e changes drastically. We mention in passing that the total reflection 
occurs at the interface between the glass sample and air, but the critical incident angle is 
determined by the refractive index of the hemisphere, as can be readily seen from Snell’s law. 
In VAE 633, we found that the critical angle is 41.3°, which is in agreement with n^^ = 1.52. 
Above 0^ =42°, both E p and E, are totally reflected and acquire different phase shift, i.e. 
|EJ « I E, | , 0 p j* 4>,. So we have e ^ 0. 


46 



IH.2 Sample PC3.5B01 and PC3.5B02 

The measured point is located at about 3 cm from the center of the disk and at the center of a 
sector with grooves. To reduce diffraction effects, the groove side of the substrate was placed 
in contact to the hemisphere and the grooves were thus filled with index matching fluid. The in- 
plane projection of the plane of incidence is perpendicular to the track direction (♦*. = 0°). 

Figures 7 (a) and (b) show the measurement results and the matched curves for sample 
PC3.5B01 with ^ = 30° and 60°, respectively. The best fit for the data set with fa = 30° 
gives n, = 1.58, (n t - n*) = 1.5 xlO' 5 , (n, - n 3 ) = 5.87X10- 4 and fa^ = -0.3°. The fitting 
error is 2.9°. The parameters found for fa = 60° are: n, = 1.58, (n, - n^ = 1.3 xlC 5 , (n,-n 3 ) 
= 5.87x10“*, and fa** = -5.3°. The average error for fa = 60° is 1.7°, smaller than that 
f° r ^ poi = 30°. We see that the two sets of fitted parameters are nearly identical. Ingnoring the 
small values of one can assert that the principal axes are in fact along the radial, track 

and perpendicular directions. 

Some of the features shown in Fig. 6 are common to all the substrates and can be easily 
explained. First, the absolute values of 6 and e are small at small 0^. The reason is that the 
beam close to normal incidence only experiences the small in-plane birefringence, but not the 
larger vertical birefringence. The second feature is that the oscillation interval, i.e. the 
difference between incident angles corresponding to two adjacent peaks, decreases with 
increasing 0^. For example, the peaks of e in (a) appear at 0^ = 22°, 37°, 49°, 60° and 68°. 
The corresponding intervals are 15°, 12°, 11° and 8°, and decrease monotonically. This feature 


could be understood as the change of optical path length per incident angle increases with 0^. 
Another feature is that 6 always shows a 90° jump when e = 45°, see the data points at 0^ = 
22° in (a). What happens here is that the polarization ellipse is very close to a circle. A slight 
decrease in the major-axis or increase in the minor-axis causes exchange of the axes and lead 
to a 90° jump. This transition is shown in Fig. 7. 

Sample PC3.5B02 is another substrate in the same series of product as PC3.5B01. We measured 
the same position (3 cm from the center of the disk and center of a sector with grooves) in both 
VAE633 and VAE780 for = 0° and ^ = 60°. The data taken with VAE633 has a shape 
similar to that shown in Fig. 6 (b). However, the best fit gives ^ m . pUne = -50.3°. In comparison 
with PC3.5B01, the in-plane birefringence is three times smaller: (n, - n^ = 0.5 xlO' 5 and 
vertical birefringence is somewhat larger (n, - n 3 ) = 6.26x1 0 4 . PC3.5B02 was also measured 
with VAE780 in the transmission mode. Figure 8 shows the measured data and the best fit. The 
best match gives = 56.7°. The fitted parameters are: n, = 1.58, (n, - n 2 ) = 1.9 xlO' 5 , 
and (n, - n 3 ) = 5.83x10^. 

m.3 Sample APO5.25B01 

The measured position is located 5 cm from the center of the disk and in the middle of a sector 
with grooves. The groove side is in contact with the hemisphere. The sample was measured with 
VAE633 at 4^ =0° for = 60°. The measured data and the best fit are shown in Fig. 9. 
Both in-plane and vertical birefringence are small: n, = 1.50, (n, - n 2 ) = 0.4 x 10' 5 , (n x - n 3 ) = 
1.3X10’ 5 and ^. pUne = -72.5°. 
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IH.4 Sample PC5.25B01 and PC5.25B02 

These two samples are from the same series of product. The measured point on smaple 
PC5.25B01 was located 5 cm from the disk center and in the middle of a sector with grooves. 
The groove side was faced to the hemisphere and filled with index matching fluid. The 
measurement was performed with fa = 0° and fa = 60° with VAE780. The measured data 
and the match are shown in Fig. 10. The fitted parameters are: n, = 1.58, (n, - = 1.3 x 10 5 , 

(n, - n 3 ) = 6.22x10“* and fa pltM = 66°. 

For Sample PC5.25B02 we measured the region without grooves: a point located 2.5 cm from 
the center of the disk with VAE633 for fa = 30°, 45°, 60° and fa = 90°. Figure 11 (a)-(c) 
show the measured data and the fits. No match has been found for 0^ >54° (not plotted). The 
fitted parameters based on data taken for 10° ^ 0^ < 53° with fa = 30°, 45° and 60° are: 

n i = 1*58, (n, - nj) = 2.2x10 ' 5 , (n ( - n 3 ) = 5.83x10^ and r< = 1.6°. We also measured 

a point in the region with grooves. This time we found that, for 0^ ^ 25°, the measured 
curves cannot be matched. The measurement was repeated for the same point, but with the side 
of grooves facing backward, so that the grooves were not filled with the index-matching fluid. 
In this case there were a series of diffracted beams and we measured the one corresponding to 
the zeroth order pattern. The measured values for 6 and e in the range of 10° ^ 0,^ < 60° 
were the same as those with the filled grooves. 

m.5 Sample PC3.5C01 and PC3.5C02 

These two MO disks are from the same series of product with quadrilayer structure. The 



measured point was located 3 cm from the disk center and in the middle of a sector with 
grooves. The data were taken with = 90° and ^ = 60°. The measured data and the best 
fit are shown in Fig. 12. The fitted parameters are: n, = 1.58, (n, - n 2 ) = 2.3 xlO' 5 , (n, - n 3 ) 
= 5.34x10'*, and = -2.9°. The principal axes are in fact along the radial and track 

directions. To obtain the best fit, the quadrilayer structure must be used and adjusted in the 
calculation. We found that the searched structure parameters were not unique. However, the 
searched birefringent parameters were always the same and hence are reliable. The reason is that 
the measured 6 and e contain information about the phase shift contributed by the birefringence, 
but not much information about the absolute reflectivities which are signatures of the 
quadrilayer. The same measurement was performed for PC3.5C02. The fitted parameters are: 
n, = 1.58, (n, - n^ = 2X1C 5 , (n, - n 3 ) = 5.51X1C 4 , and ^ hrlmf = -6.5°. Compared with the 
bare substrates PC3.5B01 and PC3.5B02, we see that the coated substrates have larger in-plane 
birefringence but smaller vertical birefringence. The decrease of vertical birefringence due to 
coating was also reported in [11], 

in.6 Sample PC5.25C01 

The sample is a small section of coated substrate with grooves cut from a 5.25" MO disk. The 
data were taken with VAE633 at = 90° for ^ = 45° and 60°. As shown in Fig. 13, both 
measured curves can be fitted very well with a fitting error of 1.65°, using the set of 
parameters: nj = 1.58, (n t - nj) = 0.9xlC 5 , and (n, - n 3 ) = 4.61x10'* and = 8°. 
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m.7 Sample PC5.25C02 and PC5.25C03 

These two samples are from the same series of product. For PC5.25C02 we measured 4 points, 
each located at the same radial position (5 cm from the disk center) but in the middle of different 
sectors. $ mc = 90° and ^ = 60° were fixed. The measured data were nearly identical for the 
4 points, see Fig. 14, indicating that the birefringence property is uniform along the track. The 
best match for one of the measured points (solid lines in Fig. 14) gives: n t = 1.58, (n, - nj = 
0.7xl0' 5 , and (n, - n 3 ) = 3.37X10 4 and ^.pu^ = -18°. 

For PC5.25C03 we measured one point at 5 cm from the center of the disk and in the middle 
of a grooved sector. The fitted parameters are: n, = 1.58, (n, - n 2 ) = 1.9 xlO' 5 , (n, - n 3 ) = 
3.47X10 -4 and \p m . pltoc = 11.5°. Figure 15 shows the measured data for both PC5.25C02 and 
PC5.25C03 with ^ = 60°, emphasizing the similarity of the two disks. 

IV. CONCLUDING REMARKS 

In this paper we have described an approach to measuring the birefringence for optical disk 
substrates using ellipsometry. The procedure consists of measurement of the rotation and 
ellipticity, and data fitting using a computer program. Using this method we have measured a 
number of bare and coated substrates. One of the main differences between our method and the 
previous ones is that we can measure the orientation of the principal axes of the index ellipsoid. 
Our results confirm that the disk normal direction is always a principal axis (e.g. tilted angle of 
less than 0.1° was found for all the measured samples). However, rotation of the in-plane axes 
was found to occur frequently. This rotation has been mentioned previously [9], but our work 



gives the first quantitative results. Another advantage is due to the use of the glass hemisphere 
which makes large propagation angles. We found that the grooves in bare (even when filled with 
index-matching fluid) and coated substrates have different effects on the polarization state of the 
beam. For bare substrates, the measured data can be matched very well when the plane of 
incidence is perpendicular to the grooves ($« = 0°). The opposite is true for coated substrates, 
i*c. , — 90 gives the best match. For other values of the match is usually poor, 

especially for data measured at large ©«. This dependence on shows that the grooves can 
cause a phase shift, thus producing errors in birefringence measurement. The fact that we can 
get good match for data taken at $>** = 0° for bare substrates and = 90° for coated 
substrates may suggest that at these particular angles the grooves have no effect on the beam, 
and thus the measured birefringence data are more accurate. 
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FIGURE CAPTIONS 

FIG.l The system setup for the variational-angle ellipsometer with wavelength X = 632.8 nm 
(The VAE633 system). 

FIG.2 (a) The system setup for the variation-angle ellipsometer with wavelength X = 780 nm 
(The VAE780 system), (b) The hemisphere setup in transmission mode, (c) The hemisphere 
setup in reflection mode. 

FIG. 3 (a) Definition of the incident polarization angle (b) The polarization ellipse of the 
reflected beam is characterized by 6^ and e. The sign of e is for LEP and " + " for REP. 
(c) and (d): Polarization state of the beam after the QWP when the fast axis is aligned with the 
major axis. 

FIG.4 (a) The quadrilayer structure of MO disk is used in fitting the data for coated substrate, 
(b) The in-plane principal axes are found rotated for some substrates, but the third one is always 
perpendicular. 

FIG. 5 The best fit for glass sample, (a) The medium is forced to be isotropic, (b) Birefringence 
is allowed. Data were taken with VAE633. The small fluctuations appearing in the calculated 
curves are due to the interference of beams reflected at the front and the back surfaces of the 
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sample. 

FIG.6 The best fit for sample PC3.5B01 measured with (a) ipp<* = 30°, (b)^ = 60°. The 
discontinuity at 0^ = 22° in (a) is explained in the text. The jumps at 0^ = 27° and 39° are 
not physical, since ±90° stand for the same polarization plane. Data were taken with VAE633. 
To obtain good match for bare substrates, we found the incident plane must be perpendicular to 
the track direction ($*. = 0°). 

FIG. 7 The exchange of major axis and minor axis causes a discontinuity of 90° in 6 when e = 
45°. 

FIG. 8 Measured data and the best fit for Sample PC3.5B02. Data were measured with VAE780 
in the transmission mode. 

FIG.9 Measured data and the best fit for Sample APO5.25B01. Data was taken with VAE633. 
The fluctuations of the calculated curves are due to the interference between beams reflected 
from the front and back surfaces of the substrate. 

Fig. 10 Measured data and the best fit for Sample PC5.25B01. Data were measured with 
VAE780 in the transmission mode. 

FIG. 11 Measured data and the best fit for Sample PC5.25B02. The measured point locates in 
a region without grooves. The polarization angle is (a) ^ = 30°, (b) ^ = 45°, (c) ^ = 
60°. No match has been found for 0^ ^ 54° (not shown). 

FIG. 12 Measured data and the best fit for Sample PC3.5C01. The data were measured at 
= 90° and ^ = 60° with VAE633. In contrast to bare substrate, we found that, for coated 
substrates, good match can only be achieved with ^ = 90°. 

FIG. 13 Measured data and the best fit for Sample PC5.25C01. The data were measured at 4*^ 



= 90° with VAE633. (a) ^ = 45°, (b) ^ = 60°. 

FIG. 14 Comparison of measured data for 4 points at different sectors of Sample PC5.25C02. 
Solid curves are for sector 1; dashed for sector 2; dotted for sector 3 and dash-dotted for sector 
4. The data were taken with VAE633 at <l' inc . = 90° and ^ = 60°. (a) Rotation angle 6, (b) 
Ellipticity s. 

FIG. 15 Comparison of measured data for the same position at PC5.25C02 (solid lines) and 
PC5.25C03 (dashed lines). The data are almost the same for 0^ ^ 45°. This shows the 
birefringent property in the two samples is very close. 
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Summary of measurement results. 
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Ring Lens Focusing and Push-Pull Tracking Scheme for Optical Disk Systems 


R. Gerber, J. Zambuto, J.K. Erwin, and M. Mansuripur 
Optical Sciences Center, University of Arizona, Tucson, Arizona 85721 


Abstract 

An experimental comparison of the ring lens and the astigmatic techniques of generating focus-error- 
signal (FES) in optical disk systems reveals that the ring lens generates a FES over two times steeper 
than that produced by the astigmat. Partly due to this large slope and, in part, because of its 
diffraction-limited behavior, the ring lens scheme exhibits superior performance characteristics. In 
particular the undesirable signal known as "feedthrough" (induced on the FES by track-crossings 
during the seek operation) is lower by a factor of six compared to that observed with the astigmatic 
method. The ring lens is easy to align and has reasonable tolerance for positioning errors. 
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Introduction : The ring lens technique of focus error detection and its combination with push-pull 
tracking was proposed in an earlier paper [1]. Later, based on diffraction calculations, the ring lens 
was predicted to have superior characteristics over other methods of focus-error detection [2, 3]. This 
superiority stems from the fact that the ring lens produces a large, yet diffraction-limited image at 
the detector, has circular symmetry which simplifies alignment procedures, and allows focus- and 
track-error signals as well as data signals to be detected in the same channel. This paper reports the 
experimental verification of some of the above predictions; in particular, we report that the ring lens 
produces a steeper FES and has significantly lower feedthrough (i.e., focus error signal generated 
when moving across tracks) than the astigmatic technique. In the sections that follow we describe a 
method that has been used successfully to fabricate the ring lens, describe our static testbed for 
making measurements of the focus- and track-error signals, present the results of these measurements 
for a ring lens, and compare the results with those obtained for the astigmatic scheme. 

Ring Lens Fabrication : The ring lens used in these experiments was formed on a soda-lime glass 
slide by a thermal process. Initially, a base plate of polished stainless steel with a circular groove 
carved on its flat surface was prepared. This plate was coated with a thin layer of black ink which 
serves as mold release agent, and then placed in contact with the glass slide. The assembly was heated 
in an oven to 735°C for 2 minutes, then slowly cooled back to the room temperature. During heating, 
the softened glass sagged into the groove and formed a curved surface in the annulus defined by the 
groove. The glass slide was subsequently removed from the mold, cleaned, and without further 
polishing was used in the experiments. 

Figure 1 shows various features of the fabricated ring lens. Dimensions of the lens are depicted 
in FigA(z) as follows : inner diameter = 2.03mm, outer diameter = 4.35mm, focused ring radius = 
1.595mm. The unshaded annular region in the figure corresponds to the section of the lens that has 
been profiled by the WYKO Corporation's TOPO system and shown in F/g.l(b). A scanning electron 
micrograph of a section of the ring lens appears in (c), and an interferogram showing Newton rings 
for this lens is reproduced in (d). 

It must be emphasized that although the thermal process described above yields reasonably good 
ring lenses, we do not consider it suitable for the fabrication of highly accurate, aberration- free 
lenses. Other shortcomings of this technique include the inability to produce lenses with very small 
ring radii, and also lack of any direct control over the focal length of the lens. Ideally, such lenses 
must be made by direct molding of glass or plastic, or by the fabrication methods of holographic 
optics or binary optics. 

Testbed for Measurements of Focus/Track Error Signals : The system used to test the various 
focusing techniques is the static focusing/tracking testbed shown in Fig. 2. This system operates at 
two HeNe wavelengths (633 nm and 544 nm), and allows either of two linear polarizations (parallel 
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and perpendicular to tracks) to be used for measurements. The paths of the red and green beams are 
coincident after their combination by the dichroic beamsplitter. The beams are then split into two 
beams of orthogonal polarizations and recombined by the polarizing beamsplitter. One can readily 
select a wavelength and a polarization direction by blocking the undesired paths on the testbed. The 
beam intensity may be controlled by a variable attenuator placed after the dichroic splitter. 

After passing through a spatial filter and a collimating lens, the beam goes through a removable 
shear plate collimation tester, which can detect wavefront errors as low as 0.1 A. The beam emerging 
from the lens is fairly large in cross-section, and is stopped by a 5 mm diameter aperture. The 
collimated beam passes through a 50/50 beam-splitter cube which sends half of the light to the disk 
(or an immitation of the disk) and the other half to a plane mirror, which acts as the reference arm in 
a Twyman-Green interferometer. The reference arm is used to check the collimation of the beam 
returning from the disk, and is covered during focus/track error measurements. 

In our experiments the grooved surface of an optical disk was simulated by a reflection 
diffraction grating, placed at the focal plane of a 40x, 0.55 NA, infinity-corrected microscope 
objective. The grating has 1.1 /im track-pitch, a sinusoidal profile, and peak-to-valley depth of 
67 nm. The choice of grating instead of an actual grooved disk eliminates the problem of focusing 
through the substrate, which would unnecessarily complicate the measurements by introducing 
aberrations due to misalignment and birefringence. The use of grating also affords us more flexibility 
in selecting the numerical aperture of the objectives, which are generally not available in a wide 
variety for 1.2 mm cover slides. The position of the grating is controlled by low-voltage piezoelectric 
transducers (PZT), mounted in the actuators of an XYZ translation stage. A GPIB -controlled power 
supply can move the PZT's in 0.003 micron increments over a range of 12 microns. 

The focusing/tracking optics are placed after the beam-splitter cube, and form an image on a 
CCD camera whose pixel spacing is 10 fim. The ring lens forms an image directly on the camera, 
while the spot formed by the astigmat is small and requires magnification. The image from the 
camera is fed into a frame-grabber card, which can simulate arbitrarily-shaped detectors in software. 
The simulated detector for the ring lens scheme is a phi detector (so named because of its 
resemblance to the Greek letter <f>), while that for the astigmatic method is a quad detector. 

Focus-Error Detection and Push-Pull Tracking with the Ring Lens : The ring lens used in our 
experiments had a focal length of 25 mm, and focused a collimated beam onto a sharp, nearly 
diffraction-limited ring of diameter 3.19 mm. A section of the focused ring is shown in Fig. 3, as an 
intensity plot of the light pattern collected by the camera; the width of this ring is roughly 50 nm. 

A schematic diagram of the ring lens focusing technique is shown in Fig. 4. As the disk (grating) 
moves in and out of focus, the beam incident on the ring lens becomes divergent or convergent, and 
the light at the focal plane of the ring lens moves outside or inside its nominal radius. A focus error 
signal may thus be generated by subtracting the integrated light intensity outside the nominal ring 


radius from that inside. The difference signal is subsequently normalized by the sum signal in order 
to eliminate the effects of laser power and/or reflectivity fluctuations. 

During seek operations when the beam scans across the tracks while remaining in focus, 
diffraction effects produce a so-called "baseball pattern" on the returning beam. This pattern 
produces bright and dark regions on the focused ring. A track-crossing signal can thus be generated 
by subtracting the integrated intensity in one half of the image from that in the other half; again the 
difference signal should be normalized by the sum for reasons that were described earlier. 

Focus-Error Detection and Push-Pull Tracking with the Astigmatic Lens : The astigmatic 

scheme is depicted in Fig. 5. The astigmat used in these experiments had focal lengths 19.75 mm and 
20.25 mm, and a clear aperture of about 5mm. The detector was placed at the circle of least 
confusion, midway between the two line foci, where deviations of the incident beam from collimation 
produced an elongated image in one of two mutually orthogonal directions. The elongation of the 
spot is detected by a four-quadrant detector; the FES and TES signals are generated by subtracting 
the integrated intensity on one pair of quadrants from that on the other, then dividing by the sum 
signal. Since in our experiments the spot itself was too small to be properly resolved by the CCD 
array, the plane midway between the foci of the astigmat was magnified by a 40x magnifier and 
imaged onto the camera. All subsequent processing of the signals was carried out on this magnified 
image. 

Results and Discussion : The ring lens technique was found to produce a steeper FES curve and 
significantly less feedthrough than the astigmatic method. These findings apply to both wavelengths 
and both polarization directions employed in this work. The ring lens was also found to be fairly 
insensitive to misalignments. 

Figures 6 (a,b) show the FES curves for the ring lens and astigmatic techniques, respectively. The 
four curves in each group are for different combinations of wavelength and polarization. Neither 
polarization nor the wavelength seem to have a significant effect on the shape and slope of the FES 
curve. The slope of the linear region of the FES for the ring lens is 0.5 per micron of defocus, 
whereas that for the astigmat is only 0.22 per micron. 

Plots of TES and spurious FES during track-crossings (i.e., feedthrough) are shown in Fig. 7 for 
the ring lens and Fig. 8 for the astigmat; important features of these plots are also summarized in 
Table I. Note that the TES amplitudes are comparable for the two schemes. On the other hand, for 
both wavelengths and polarizations, the feedthrough (in microns) for the ring lens is between six and 
eight times lower than that for the astigmat. 

The ring lens technique is fairly insensitive to misalignments and imperfections of the system. 
Figure 9(a) shows several FES curves for various radii of the phi detector. (Recall that the radius of 
the focused ring is around 1.6 mm, and its width is approximately 50 pm.) An error of ± 10 fim in 
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detector radius will change the zero-crossing by ± 1.6 microns. In an actual drive, if the radius of the 
detector happens to deviate from that of the focused ring, the electronic gains of the detectors must 
be adjusted to offset this error. Figure 9(a) also shows that a mismatch of ± 20 pm between the ring 
radius and the radius of detector can seriously degrade the FES slope. In Fig. 9(b) we show FES 
curves for the best radius but various locations of the phi detector in the focal plane of the ring lens. 
These plots show a slight decrease in the slope as well as a small shift of the zero-crossing for 10 nm 
of centering error. The consequences of a 20 nm decenter, however, are much more drastic. 


Focusing Wavelength Polarization Peak-to-valley TES Peak-to-valley feedthrough 


Scheme 

(nm) 

(relative to tracks) 

norm, units 

norm, units 

microns 

wavelengths 

Ring lens 

544 

parallel 

1.9 

0.12 

0.25 

0.46 

Ring lens 

544 

perpendicular 

1.9 

0.12 

0.25 

0.46 

Ring lens 

633 

parallel 

1.5 

0.07 

0.14 

0.23 

Ring lens 

633 

perpendicular 

1.5 

0.08 

0.17 

0.27 

Astigmat 

544 

parallel 

1.6 

0.34 

1.5 

2.75 

Astigmat 

544 

perpendicular 

1.7 

0.44 

2.0 

3.65 

Astigmat 

633 

parallel 

1.4 

0.22 

1.0 

1.60 

Astigmat 

633 

perpendicular 

1.4 

0.22 

1.0 

1.60 


Table I. Summary of track-crossing and feedthrough signals 

In comparison with the ring lens results in Fig. 9(b ) , Fig. 10 shows several FES curves for various 
locations of the quad detector of the astigmatic technique. We observe that 6 nm of centering error 
has negligible effects on the FES curve, but a decenter of 12 fim is beginning to cause severe 
problems. 

Concluding Remarks : We have presented a comparison of the ring lens and astigmatic methods of 
generating focus- and track-error signals for optical disk drives. It was found that the ring lens 
produces an FES curve which is substantially steeper than that produced by the astigmat. Also the 
cross-talk on the focus error signal from track-crossings was 6-8 times smaller for the ring lens 
method. In addition, we found the ring lens scheme to be tolerant of small errors in detector 
positioning and alignment. 
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Figure Captions 




Fig. 1. a) Diagram showing dimensions of the ring lens. The apex of the lens occurs at a radius of 
1.595mm, which corresponds to the focused ring radius of a collimated incident beam. 

b) Surface profile of section of the lens obtained by interferometric measurements. The 
region captured by the profilometer is the unshaded annular region depicted in (a) with a 
width of 0.524mm. 

c) SEM image of section of the ring lens. 

d) Newton's rings observed on a section of the ring lens. 

Fig. 2. Schematic diagram of the static focusing/ tracking testbed. 

Fig. 3. Section of the intensity profile of the focused ring (pixel = 10 urn x 10 urn). 

Fig. 4. Schematic diagram of the method of ring lens focusing and push-pull tracking. 

Fig. 5. Schematic diagram of the method of astigmatic focusing and push-pull tracking. 

Fig. 6. Measured FES curves for the ring lens (a) and astigmatic (b) scheme of focus error 
detection. The various symbols in each plot correspond to different wavelengths and 
polarization orientations used in the experiment. 

\ 

Fig. 7. Measured TES (dashed) and feedthrough (solid) curves for the ring lens scheme, 

a) Green laser with polarization parallel to tracks, 

b) Green laser with polarization perpendicular to tracks, 

c) Red laser with polarization parallel to tracks, 

d) Red laser with polarization perpendicular to tracks. 

Fig. 8. Measured TES (dashed) and feedthrough (solid) curves for the astigmatic scheme. 

a) Green laser with polarization parallel to tracks, 

b) Green laser with polarization perpendicular to tracks, 

c) Red laser with polarization parallel to tracks, 

d) Red laser with polarization perpendicular to tracks. 
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Fig. 9. FES curves for the ring lens scheme in the case of red laser with polarization parallel to 
tracks. 

a) Different curves correspond to different detector radii. The solid curve is obtained for 
best radius, V is for ± 10 nm and "x" for ± 20 pm change in the radius. The curves shift to 
the left with increasing detector radius. 

b) Different curves correspond to various detector locations within the focal plane of the 
ring lens. The solid curve ("o") is for the best detector position, while the dashed curves are 
for translations of 10 \im (V) and 20 \im ("x") relative to this optimum position. 

Fig. 10. FES curves for the astigmatic scheme in the case of red laser with polarization parallel to 
tracks; different curves correspond to various detector locations within the midplane 
between the two astigmatic line foci. The solid curve ("o") is for the best detector position, 
while the dashed curves are for translations of 6 pm ("+") and 12 pm ("x") relative to this 
optimum position. 
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Figure 3 
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FES curves for astigmat, both wavelengths and polarizations 
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Tracking offset in microns 
Figure 7(a) 
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ABSTRACT 


The domain structure of Co/Pd (1.6/6. 3 A)xN multilayers and its relation to the bulk 
magnetic properties of the samples were studied. The Co/Pd multilayeres were deposited by rf 
and dc magnetron sputtering onto different substrates. It was found that magnetic and magneto- 
optical properties and domain structure of the multilayers were affected by total film thickness 
and substrate condition. Magnetization, coercivity and anisotropy of the films decreased 
significantly as the film thickness dropped below 100 A. However, Kerr rotation angle had a 
maximum at the same thickness. The width of the domain structure increased with the decrease 
of the film thickness attaining the single domain state at N=10. The initial curves in Co/Pd 
multilayers were found to depend on demagnetization process. The samples demagnetized by in- 
plane field showed the largest difference between initial curves and the corresponding parts of 
the loops. Different domain structures were observed in the samples demagnetized by 
perpendicular and in-plane magnetic fields. 
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I. INTRODUCTION 


Multilayered films of Co/Pt and Co/Pd are promising materials for magneto-optical (MO) 
recording [1,2]. Compared to amorphous rare earth-transition metal alloys [3] currently used in 
commercial MO disks the multilayers are attractive for their superior corrosion resistance [4] 
and larger Kerr effect signal at short wavelengths [5] (<500 nm), where higher density 
recording is possible. Co/Pt multilayers have generally been considered more attractive than 
Co/Pd because of their larger Kerr rotation angle. Co/Pt is deficient, however, in coercivity. 
The coercivity of Co/Pd multilayers is considerably higher than that for Co/Pt [6]. The 
coercivity plays an important role in the recording process, since it governs the domain size, 
shape and stability during readout. Any fluctuations in coercivity may cause an inhomogeneity 
in domain shape, which in turn contributes to recording noise. The study of coercivity 
mechanisms in the multilayered structures is also of great importance from the fundamental point 
of view. The coercivity is closely related to magnetic anisotropy of thin films as well as to the 
microstructure parameters of the magnetic media such as interface quality, growth orientation, 
crystallite size, distribution of easy axis for individual crystallites and so forth. The shape of 
magnetic hysteresis loops, initial curves and domain structure also strongly depend on 
microstructure of the films. Consequently, the investigation of domain structures in multilayered 
films and their relation to hysteresis loops and initial curves can help in understanding the 
origins of coercivity in these films. 

The present work is devoted to investigation of domain structures in Co/Pd(l .6/6.3 A)xN 
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multilayers of varying thickness fabricated on different substrates, and their relation with 
hysteresis loops and initial magnetization curves. 


H. EXPERIMENTAL 


Multilayers of Co/Pd were fabricated by sputtering in Ar at a pressure of 20 mTorr onto 
various substrate conditions. The Co and Pd layers were deposited by rf and dc magnetron 
sputtering, respectively. Four different conditions were chosen, namely the multilayers were 
deposited onto Si(lll) substrate, Si(lll) etched (80A) by Ar, unetched SiN film which was 
deposited onto Si(lll) and etched SiN film on Si(lll). The thickness of SiN film was 850 A. 
For all the cases, individual layer thicknesses of Co and Pd were 1.6 A and 6.3 A, respectively. 
During deposition the substrate was kept at ambient temperature. The Pd layer was the first 
deposited onto substrate. The number of Co/Pd pairs, N, was 50 for all types of the substrates, 
and 6, 10, 20, 30, 40 for the etched SiN substrate (see table 1). Last letters in the sample ID 
mean substrate condition, and a number before the letters is attributable to the number of Co/Pd 
pairs. For example, sample ID S30ESiN means that substrate is etched SiN layer and the film 
has 30 Co/Pd pairs. 

The measurement of magnetic properties of the samples was carried out by vibrating 
sample magnetometry in fields up to 14 kOe. The area of the samples for calculation of 
magnetization was measured using an optical microscope with digitizing image capability. 
Magneto-optic loops and initial curves were measured in our magneto-optic loop tracer [7] in 
fields up to 20 kOe. The constants of perpendicular anisotropy were determined using data 



obtained from the MO loop tracer with in-plane applied field [8]. Domain structures were 
observed with the help of polarizing optical microscope using a lOOx objective (oil immersion 
type). The observations were recorded with a PC-based frame grabber that allowed image 
processing for noise reduction, image segmentation and image thresholding [9]. A single-pole, 
conical tip electromagnet capable of producing a maximum field of 5 kOe, was mounted under 
the microscope’s XY stage to provide the necessary fields for domain growth and contraction. 

m. EXPERIMENTAL RESULTS AND DISCUSSION 
A. Domain structure and substrate condition 

Our observations of domain structure in Co/Pd multilayers indicate that shape of the 
domains and their growth strongly depends on the substrate conditions and the total thickness 
of the film. In Fig.l are shown the domain structures developing in (Co/Pd)x50 multilayered 
films fabricated on different substrates. The films S50ESi, S50USiN and S50ESiN show domain 
wall motion dominated magnetization process. Domain structure in the films initially saturated 
by the field of 20 kOe starts to develop from randomly distributed points. These points are 
believed to be some defects which facilitate creating a domain wall. Expansion of the domain 
in the magnetic field close to the coercive field looks like development of maze-like structure 
with very narrow domains. This behavior proves that DW coercivity is smaller than nucleation 
coercivity. The domain structure in the film deposited onto unetched Si substrate is quite 
different (Fig. Id). The remagnetization in this film occurs due to nucleation process. With the 
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increase of the magnetic field the entire film is covered with domains of opposite magnetization; 
thus the wall coercivity in this film is greater than nucleation coercivity. The magnetic properties 
of (Co/Pd)x50 multilayers are shown in table 1. It follows from this table that the 
macromagnetic parameters like M, and H c are very similar. K„ is about two times smaller for 
samples fabricated on unetched substrates. The lower value of K„ in the films on unetched 
substrates allows to demagnetize the samples with in-plane field of 20 kOe. The same in-plane 
field, applied to (Co/Pd)x50 films, fabricated on etched Si and SiN substrates, is not able to 
demagnetize the samples. 

The authors of [10] measured the roughness of etched and unetched Si and SiN substrates 
as well as microstructure of Co/Pt multilayers fabricated on the substrates. They showed that 
the dispersion angle for < 1 1 1 > orientation of crystallites has the largest value for samples on 
unetched Si substrates. The presented in [10] magneto-optic loops indicate that the magnetization 
process in Co/Pt film on unetched Si substrate is nucleation dominated as in our case for Co/Pd 
film on the similar substrate. So, we can assume, that general microstructural properties of 
Co/Pt and Co/Pd multilayers on Si and SiN substrates are similar. Hence, the microstructure of 
the samples seems to be responsible for the type of magnetization process. The film with low 
perpendicular anisotropy and big dispersion angle of easy axes shows nucleation dominated 
magnetization. Samples with high anisotropy and low dispersion angle of easy axes show wall- 
motion dominated magnetization process. 


B. Domain structure and film thickness 



Figure 2 shows the magnetic and magneto-optic properties of (Co/Pd)xN films of 
different thicknesses, deposited onto etched SiN substrates, as functions of film thickness. The 
magnetization of the samples remains the same for N from 50 to 20 decreasing slightly for 
N=10 and 6. The coercivity decreases starting from 40 layers and reaches very low value of 
85 Oe for N=6. Kerr rotation has a maximum for N=10, which is attributable to interference 
effects [11]. Anisotropy constant decreases with the decrease of film thickness, approaching 
a value of 7x10 s erg/cc for N=6. Especially big drop in K„ value occurs at film thicknesses 
below 120 A. The authors of [11] showed that this region of film thicknesses corresponds to 
formation of continuous film structure. For a Co/Pt film with the thickness t of 30 A the distinct 
island structure was observed. On the contrary, the film with t= 150 A had a continuous 
structure which was formed of very fine crystalline grains with the grain size below 100 A in 
diameter. 

Figure 3 represents the domain structures in (Co/Pd)xN multilayers of different 
thicknesses deposited onto etched SiN sublayer. With the decrease of the total thickness of the 
films the domain width of maze-like structure increases approaching the single domain state at 
N=10 and dropping again for N=6. The dependence of domain width on film thickness is 
shown in Fig.4. It is known from the domain theory [12] that the domain structure of the films 
with uniaxial magnetic anisotropy is effected by film thickness, t, and, in very thin films, single 
domain structure is attained. In a uniform film with strong perpendicular anisotropy the critical 
thickness at which the transition from a stripe domain to a single domain occurs is estimated by 
the formula [11]: 

t,= 16.8 <t/t 4 M, 2 
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where a=4(AK u ) 05 is domain wall energy density, M, is saturation magnetization, A - exchange 
stiffness, K„ - anisotropy constant. Using experimental data for K„ and M, and assuming A = 10' 7 
erg/cm we can estimate t c at about 200 A. Experimental results indicate that in Co/Pd 
multilayers this value is about 100 A. This discrepancy may be attributed to the fact that in the 
multilayered structures the effective exchange stiffness constant A can be lower than in pure 
ferromagnet because of the inclusion of nonmagnetic layers of Pd. The appearance of maze-like 
structure again in (Co/Pd)x6 multilayers can be explained by significant decrease in constant 
(Fig.2a), so t. is becoming lower than actual thickness of multilayered film with N=6. 

C. Domain structure and initial curves 

Measurements of the magneto-optic loops and initial curves have shown that the shape 
of the initial curve in Co/Pd multilayers strongly depends on the demagnetization procedure 
(Fig. 5). If we demagnetize the sample by a perpendicular field of the same direction as for 
measuring the initial curve, the initial curve quickly approaches the corresponding point on the 
magneto-optic loop and goes along it up to saturation value. If we demagnetize the sample by 
a perpendicular field of opposite sign or by in-plane field the initial curve goes faster than the 
corresponding part of the magneto-optic loop and approaches the saturation value in lower 
magnetic fields. This difference between loop and initial curve is greater for samples 
demagnetized by in-plane magnetic field. In Fig. 6 are shown domain structures in SlOESiN and 
S06ESiN films developed both by perpendicular and in-plane magnetic fields. Usually if we 
demagnetize by a perpendicular field, the width of the domains is comparable to that for domain 


structure expanding in the same sample . However if we demagnetize by in-plane fields, the 
width of domains is noticeably smaller than that for developing domain structure. The most 
striking difference is observed for SlOESiN film (Fig.6a,b) where single domain state is 
attained. After demagnetization by in-plane field of 20 kOe the film is covered by submicron 
domains. 

As mentioned earlier, the initial curves for samples demagnetized by perpendicular fields 
of opposite directions are somewhat different. To determine the source of this difference we 
performed the following experiment: First, the domain structure was allowed to develop in the 
magnetic field in one direction (see Fig. 7a). Then the reverse magnetic field was applied for 
several seconds to contract the structure (Fig.Tb). The values of both magnetic fields were fixed. 
Finally, the difference of the two images was obtained ( Fig.7c). The experiment has shown that 
the required field to start shrinking the structure in Fig. 7a is lower than the opposite field needed 
to expand the structure. Figure 7c shows that the contraction occurs along the length of the 
walls, rather than at domain tips. 

The difference in the initial curves for different demagnetized states can originate from 
pinning of DW on inhomogeneities in the films [13]. If we continue to expand a domain wall 
we need more field the DW to overcome energy barriers caused by a chain of pinning sites, 
which the DW has already approached. If the domain shrinks, the DW has some freedom to 
move between adjacent chains of pinning sites. So, the field required to shrink the DW is lower 
than that to expand the domain. Demagnetization by in-plane field produces much smaller 
domains. They arise randomly without any correspondence to pinning sites. In this case the 
domain walls have even more freedom to move. Moreover, small diameter facilitates collapse 


116 


of the domains. Therefore the initial curve for the films demagnetized by in-plane field goes 
much faster than that in the films, demagnetized by perpendicular field. 


IV. CONCLUSION 


It was found that magnetic and magneto-optical properties and domain structure of Co/Pd 
multilayers were affected by total film thickness and substrate conditions. Both wall motion- and 
nucleation-dominated processes were observed in the samples fabricated on different substrates. 
Magnetization, coercivity and anisotropy of the films decreased significantly as the film thickness 
dropped below 120 A. The width of the domain structure increased with the decrease of film 
thickness, attaining the single domain state at N=10. The initial curves in Co/Pd multilayers 
were found to depend on the process by which demagnetization had been attained. The samples 
demagnetized by in-plane field showed the largest difference between initial curve and the 
corresponding parts of the loop. Different domain structures were observed in the samples 
demagnetized by perpendicular and in-plane magnetic fields. 
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Table 1: Co/Pd(1.6/6.3 A)xN samples 
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Figure la. Domain structure in Co/Pdx50 multilayered film, 
deposited onto etched SiN layer, developing at H=2.8 kOe. 
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Figure lb. Domain structure in Co/Pdx50 multilayered film, 

deposited onto unetched SiN layer, developing at H=2.29 kOe. 
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Figure lc. Domain structure in Co/Pdx50 multilayered film, 

deposited onto etched Si substrate, developing at H=2.8 kOe. 
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Figure Id. Domain structure in Co/Pdx50 multilayered Film, 

deposited onto unetched Si substrate, developing at H=2.67 kOe. 
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Figure 2b. Dependence of Kerr angle and coercivity on the 
thickness of Co/Pd multilayers. 
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Figure 3a. Domain structure in SSOESiN sample developing at H=2.8 kOe. 
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Figure 3b. Domain structure in S40ESiN sample developing at H=2.8 kOe. 
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Figure 3c. Domain structure in S30ESiN sample developing at H=2.44 kOe. 


l- 10/zinH 



Figure 3d. Domain structure in S20ESiN sample developing at H=1.63 kOe. 
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Figure 3e. Domain structure in SlOESiN sample developing at H=0.85 kOe. 
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Figure 4. Dependence of average domain width on the thickness 
of Co/Pd multilayers. 
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Figure 5a. Sample SlOESiN. Kerr rotation hysteresis loop and initial curve with 
perpendicular demagnetized field Hd=+580 Oe. 
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Figure 5b. Sample SlOESiN. Kerr rotation hysteresis loop and initial curve with 
perpendicular demagnetized field Hd=-620 Oe. 
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Figure 5c. Sample SlOESiN. Kerr rotation hysteresis loop and initial curve with 
in-plane demagnetized field. 





Figure 6b. Domain structure in SlOESiN sample after demagnetization by in-plane 
magnetic field Hd=20 kOe. 
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Figure 6c. Domain structure in S06ESiN sample developing at H=0.12 kOe. 

10/iITH 



Figure 6d. Domain structure in S06£SiN sample after demagnetization 


by in-plane magnetic field Hd=20 kOe. 
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Figure 7a. Domain structure in S20ESiN sample developing at H=1.63 kOe. 
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Figure 7b. Domain structure in S20ESiN sample after application of reverse field. The 


initial domain structure is shown in figure 7a. 
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Determination of Average Refractive Index of Spin Coated DCG Films 

for HOE Fabrication 

T. J. Kim, Eugene W. Campbell, and Raymond K. Kostuk 
Department of Electrical and Computer Engineering 
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Tucson, AZ, 85721 


Abstract 

The refractive index of holographic emulsions is an important parameter needed for design- 
ing holographic optical elements (HOEs). In this paper theoretical calculations of the accuracy 
required for the refractive index and thickness of emulsions needed to meet predetermined Bragg 
angle conditions are presented. A modified interferometric method is used to find average refrac- 
tive index of the unexposed and the developed dichromated gelatin holographic films. Slanted 
transmission HOEs are designed considering the index and thickness variations, and used to verify 
the index measurement results. The Brewster angle method is used to measure surface index of 
the unexposed and the developed films. The differences between average index and surface index 
are discussed. Theoretical calculation of the effects of index variation on diffraction efficiency, and 
experimental results for index modulation variation caused by process changes are also presented. 


l.Q Introduction 

The refractive index of dichromated gelatin (DCG) film is a function of many parameters. These 
include gelatin type, dichromate concentration, film thickness, environmental conditions (relative hu- 
midity and temperature), hologram construction conditions (exposure energy, construction geometry, 
wavelength) and chemical development. A number of researchers have previously investigated the 
refractive index of DCG films[l]-[8j. Several different techniques were used to investigate the index, 
such as Abbe refractometry, Brewster angle methods, interferometry with index matching oil, and 
diffraction efficiency measurements with coupled wave analysis. The reported index values range from 
1.20 to 1.71. This large uncertainty in index makes it difficult to accurately design holographic optical 
elements. 

In order to design HOEs properly, the average refractive index before exposure and after devel- 
opment must be known within an acceptable range. This paper presents the experimental results for 
the surface and the average refractive index of spin coated DCG holographic films. It also shows how 
variations in these values affect the reconstruction angle for slanted volume holograms. 


2.0 The Effect of Index Change on Bragg Angle 

DCG film can swell or shrink after chemical development. This causes the grating K vector to 
rotate. The Bragg reconstruction angle (or angle for maximum diffraction efficiency) caused by a shift 
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in the K vector can be predicted using coupled wave analysis[9]. If the emulsion thickness during 
construction is t c and the thickness during reconstruction is t r a thickness factor can be defined as 


m — tf^t^* 


The K vector will change according to 


I< r =K C . 

— r* — c x 


(i) 


( 2 ) 


where the subscripts r and c denote the reconstruction and construction, respectively. The slant angle 
cf> will be modified as 


tand> r = m • tan<£ c . 


(3) 


From Bragg angle condition, cos (<f> — y?) = | K | /2/?[9], the reconstruction angle in the medium can be 
determined from 


<^ r = d>r + COS - 7T, (4) 

where /3 r = 27rn r /A r . The shift in the Bragg reconstruction angle is a function of the thickness change 
(m), and the average refractive index (n r ) after development. It varies with the value for the K vector 
slant angle as shown in Fig. 2 (10°), Fig. 3 (30°), and Fig. 4 (50°). Indices of the surface and the 
bulk area for the unexposed film were almost the same which were measured using the Brewster angle 
method (see 3.1) for the surface and using a modified interferometric method (see 3.2) for the bulk 
area. 

When the K vector slant angle increases, the shift in Bragg reconstruction angle becomes larger in 
air for the same average refractive index variation. For example, when the desired reconstruction angle 
is 9.3° ± 0.5°, the predicted average refractive index of n r = 1.44, n c = 1.59, and the thickness swells 
by 5 % (m=1.05), the allowed average refractive index can range from 1.27 to 1.60 (Fig.2). However, 
when the desired reconstruction angle is 45.4° ± 0.5° with n r = 1.44, n c = 1.59, and m=1.05, the 
allowed average refractive index range is reduced to 1.415 to 1.465 (Fig.4). Therefore, the range in 
average refractive index must be reduced for gratings with higher K vector slant angles. 

3.0 Refractive Index Measurement Techniques 

All the DCG spin coated films used for these experiments have a composition consisting of [Ammo- 
nium dichromate : Gelatin : DI water = 10 gram : 30 gram : 250 ml]. The film thickness was measured 
with a Tancor alpha-step 200 with a vertical resolution 0.005 /tm. The unexposed DCG film is isotropic 
and uniform (t c ± 0.12 /rm) over 10 mm. Thus, it is relatively easy to measure the refractive index. 
However, it is very difficult to measure the refractive index for the developed hologram because the film 
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surface is not as uniform (t r ± 0.32 /im) and the phase variation through the grating is not linear[7], 
[8]. Therefore, several refractive index measuring techniques were used to evaluate developed DCG 
holographic gratings. The Brewster angle measurement technique and Fabry Perot resonant method 
are used to evaluate for the surface index, a modified interferometric method is used to determine the 
bulk average index of nongrating region, and the average index of grating region is measured using an 
analysis of the resultant grating efficiency. 

3.1 Brewster Angle Method 

This method (Fig. 5) was developed by Pawluczyk et al. [4], [7], [8]. It gives only the surface refrac- 
tive index[10]. Special care is required in measuring the Brewster angle because the interference of 
light reflected from air-emulsion and emulsion-substrate exists and it becomes dominant around the 
Brewster angle[4], [8]. 

The measurement results are, n c (A = 514. 5nm) = 1.590 ± 0.005, n c (A = 632. 8nm) = 1.570 ± 0.005 
for the unexposed film and n r (A = 514. 5nm) = 1.570 ± 0.005, n r (A = 632. 8nm) = 1.540 ± 0.005 for 
the uniformly exposed and developed film. 

3.2 Interference Method 


This method is very useful for measuring the bulk index of the emulsion, however, the cutting ruins 
the film (Fig.6). Calibrated Cargille index matching oil is used and the remaining index matching oil 
between trials is removed with ethyl ether. There is some degree of error if the oil is not completely 
removed between measurements. 

The film is rotated during the measurement by an angle 0 e (Fig. 6) in order to eliminate ambiguities 
in fringe matching. If the rotating angle inside the medium is 0 e and the index difference between 
emulsion and index matching oil is An, then the relative fringe movement as a function of 0 e follows 
the relationship 


2^ Ant 
Acos0 e 


= 2?r(q + 1), 


(5) 


where t is the film thickness, q is integer, and I is the fraction of relative fringe shift. When the 
shift completes one cycle by increasing the angle 0 e , I = 1, and it is possible to keep track of | An | 
variation. When the index of applied matching oil is higher than that of the emulsion, the fringes 
will either shift up or down with increasing Q t . When the applied matching oil index is lower than 
that of the emulsion, the fringe shifts in the opposite direction as 0 t increases. However, Eq. (5) can 
not be applied in the holographic grating region since the phase variation through grating does not 
vary linearly with changing 6 e . This method needs more improvements for the index measurement of 
gratings. 


Measured values for the unexposed film were n c (A = 632.8nm) = 1.58 ± 0.01, for the unexposed 
and developed film areas n r (A = 632. Snm) = 1.46 ± 0.01, and n r (A = 632. 8nm) = 1.48 ± 0.02 for the 
uniformly exposed and developed film areas. 


3.3 Grating Efficiency Analysis 

In this approach the shift in Bragg angle of a volume grating is compared to data using coupled 
wave analysis [9]. In previous work[6] polarization selective gratings were used to determine the index 
n r , and was found to be in a range from 1.2 to 1.4. The specific value depends on the emulsion thickness 
and processing technique. This experiment was repeated in our lab with an interbeam construction 
angle of 2 0 o = 30° within the emulsion. This results gave an index n r = 1.38. 

By measuring the Bragg angle shift and thickness change factor (m), the average refractive index 
can be deduced from the angular tolerance calculation as discussed in section 2.0. Transmission grat- 
ings with large slant angles (50° in air) were used since they produce larger shifts in Bragg angle. Pre 
and post processed film thickness measurements were obtained from different 2 x 2.5 inch samples 
taken from the same 4x5 inch 2 plate. The plates were dehydrated with isopropanol/water concen- 
trations of (25-50-75-100-100 %). The maximum efficiency of the hologram used for Bragg angle shift 
measurements was 8S.4 %, The Bragg angle shift was determined from the change in the center of the 
FWHM values of the efficiency vs. angle measurements and found to change by 5° (see the Exp. data 
curve in Fig. 7). This corresponds to an average index n r = 1.360 (Fig. 4). 

Theoretical calculation based on the rigorous coupled wave analysisfll] was performed and com- 
pared with the experimental data in Fig. 7. The measured thickness (t r ) is 8.698 fj.m, however, the 
effective thickness for the curve fit is 5.0 fim. This lower effective thickness may be caused by the grad- 
ual change of index modulation with depth which may be caused by nonuniform chemical dehydration. 
The difference of index between surface (1.54) and the average value (1.46) is a good evidence for this. 
More research needs to be done to explain this effect. 

3.4 Fabry Perot Resonant Method 

Another method examined to determine the refractive index of the emulsion was by comparing 
the measured reflected beam intensity to calculated reflection values. In the case of a uniform thin 
emulsion on glass the interference of beam reflected from the air-emulsion and emulsion-glass surfaces 
give rise to a Fabry Pero effect[12]. Information from the modulated beam intensity can be used to 
determine the effective index of the emulsion. However, the DCG samples used in our experiment had 
a high level of surface reticulation in the exposed areas. This resulted in high levels of scatter from 
the air-emulsion surface which obscured the Fabry Perot effect. 

In spite of these difficulties the average index of the emulsion surface could be determined, and was 
found to be 1.54 for a uniformly exposed region of film and 1.50 for a grating region at 632.8 nm. 


4.0 Effect of Dehydration Rate on Index Modulation 

The index modulation depends on the exposed energy and the chemical process. Gradual dehydra- 
tion with isopropyl alcohol concentration of (25-50-75- 1 00- 1 00%) were compared to a rapid dehydration 
(50-100-100%). 
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The 1st peak of diffraction efficiency (DE) occurs with an exposure energy of about 70 mJ/cm 2 for 
the rapid dehydration and with approximately 220 mJ/cm 2 for the gradual dehydration. Normalized 
curves for DE vs. angle are given in Fig. 8. All measurements were performed as explained in 3.3. For 
rapid dehydration the Bragg angle shift is about 7.08° and the thickness factor is about m = 1.07. The 
resultant average refractive index n r is about 1.348. For the gradual dehydration process the Bragg 
angle shift is about 5.10°, and the thickness factor is about m = 0.992. This gives an average refractive 
index n r of about 1.355. The angular band width at full width - half maximum is A0 — 13° for the 
gradual dehydration process and A 6 = 9° for the rapid dehydration. This implies that the effective 
thickness of the grating is greater for the rapidly dehydrated emulsion. 

5.0 Conclusions 

In this report several techniques for measuring the change in the average refractive index were 
evaluated. The effect of inaccuracies in values used for the pre and post processed average index on 
the diffraction angle from transmission gratings were also determined. 

The grating efficiency analysis method was the only method found which give consistent results for 
the effective index of holographic gratings. The surface refractive index for the holographic gratings 
was found to be considerably higher than the effective index (i.e. 1.50 vs 1.35). The effect of changing 
the rate dehydration did not appear to significantly change the effective refractive index. 
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Figure 1: Schematic diagram for the K vector shift where r denotes reconstruction and c construction. 
A r and A c are grating spacings. 



Figure 2: The Bragg angle shift with different thickness change factors as a function of average refractive 
index. Values given are for a slanted transmission grating; 0 = 0°, —<p = 10° in air, A c = A r = 514. 5nm. 
The point • denotes n c = n r = 1.59 and no thickness variation. 
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Figure 3: The Bragg angle shift with different thickness change factors as a function of average refractive 
index. Values given are for a slanted transmission grating; 9 = 0 °, —<p = 30° in air, A c = A r = 514. 5nm. 
The point • denotes n c = n r = 1.59 and no thickness variation. 



Average refractive index 


Figure 4: The Bragg angle shift with different thickness change factors as a function of average refractive 
index. Values given are for a slanted transmission grating; 9 = 0°, —<p = 50° in air, A c = A r = 514. 5nm. 
The point • denotes n c = n r = 1.59 and no thickness variation. 
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Longer than 1 m to reduce angular misreading 



Figure 5: Optical set-up for Brewster angle measurement technique. 



Figure 6: Optical set-up of modified interferometric technique for index matching measurement. 






Figure 7: Comparison between experimental data and theoretical calculation. Values given are for a 
slanted transmission grating; 0 = 0°, — <p = 50° in air, m = 0.998; construction and reconstruction at 
A c = A r = 514. 5nm. 



Figure 8: Experimental data of angular bandwidth for the comparison between rapid dehydration 
(Isopropyl alcohol concentration 50-100-100 %) and gradual dehydration (25-50-75-100-100 %). The 
exposure energy to obtain maximum diffraction efficiency is about 70mJ/cm 2 for the rapid dehydration 
process and about 230 mJ/cm 2 for the gradual dehydration process. Values given are for a slanted 
transmission grating; 6 = 0°, —<p = 50° in air; construction and reconstruction at A c = A r = 514. 5nm. 
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Abstract 

A shading band can alter the system transfer function. In our theory and experiment, both contast 
improvement and enhanced carrier-to-noise ratio are investigated. 
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1. Introduction 

Efforts to improve the density of data recorded onto optical media have included shaping the 
transfer function of the optical system. Several authors have described various methods for shaping the 
transfer function. Reference 1 contains a summary of the various techniques and plots of the transfer 
functions. In this paper, we describe the effects of a simple shading band in the data path of an optical 
drive. One effect of the shading band is to shape the transfer function of the optical system. This is 
similar to the technique of applying electronic equalizing filters for better frequency response. A second 
effect of the shading band is to decrease the noise associated with the laser and the medium. In our work, 
we implement a shading band in a commercial optical disk drive. The relevant parameters to this study 
are: shading band width = 30% of the pupil diameter, medium velocity = 5.7 m/s, NA = 0.55, and X = 
0.78 nm. We used a standard 130 mm magneto-optic disk with a plastic substrate. The following sections 
briefly describe the theory associated with using the shading band and the results of our experiments. 

2. Theory 

The data pattern can be considered a collection of gratings. Large bit spacings correspond to low 
frequencies and small bit spacings correspond to high frequencies. When the laser beam is focused onto 
the medium, the reflected light contains diffracted orders from the data-pattern gratings. The signal light 
is contained primarily in the ± 1st diffracted orders reflected from the disk. The amount of signal power 
that is retrieved by the data detectors depends on how much signal light is passed through the pupil of the 
objective lens. With reference to Figure 1, at low (2 Mhz) frequency the ± 1st orders are almost 
completely contained within the pupil, so signal amplitude is high. The shading band reduces the signal 
light by reducing the amount of ± 1st orders passed by the pupiL As the recording frequency is increased, 
thereby reducing the bit spacing, the ± 1st orders become more widely separated and the shading band 
mainly reduces the background light, and hence reduces the noise, in the data signal. 

The calculated transfer functions for our system with and 
without a shading band are shown in Figure 2. In our calculations, we 
assume that the data pattern consists of wide marks in the direction 
perpendicular to the tracks. When the shading band is used, the signal 
amplitude is reduced at low frequencies, but not at high frequencies. 

This modification of the transfer function partially equalizes the 
frequency response of the system, thereby providing improved contrast. 

Reduction of the signal is detrimental unless the noise is reduced as 
well. Our experiments, which are described in the next section, were 
designed to investigate the properties of signal and noise with the 
shading band. 

3. Experiment 

The optical drive was modified to accept the shading band in 
its data path. The most convenient location of the shading band is 
after the partially polarizing beam splitter and before the lens that 
focuses light onto the data detectors. 

Figure 3 displays the noise behavior of the system with and 
without the shading band. The electronic noise floor was measured 



-1 order 


+1 order 



2 MHz 


«; AHz 


6 MHz 


Figure 1. Pupil of the objective 
lens, shading band, and ± 1st 
diffracted orders of the data 
signal for three recording 
frequencies. 
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with the data detectors blocked. It is not a 
function of the shading band. The detectors 
were unblocked, and the disk was stopped to 
measure the laser noise floor. Note that the 
laser noise includes the electronic noise 
contribution. With the shading band, the 
laser noise is only 1 dB above the electronic 
noise floor. Then, with the disk spinning, we 
measured the medium noise floor on an 
erased track. The medium noise is reduced 4 
dB at low (0 to 1 MHz) frequencies, and 1 to 
2 dB at high (5 to 7 MHz) frequencies. 

Figure 4 displays the narrow-band 
carrier and noise behavior. In the low (0 to 3 
MHz) frequency region, the carrier is reduced 
about 2 dB more than the noise, resulting in a 
decrease in the carrier-to-noise ratio (CNR). 

In the mid (3 to 5 MHz) frequency range, the 
carrier and noise are reduced equally with no 

loss in CNR. In the high (5 to 6 MHz) frequency range, the carrier is not reduced, the noise is reduced 
significantly, and the resulting CNR is increased about 2 dB. 


200 


400 600 800 1000 1200 1400 1900 1800 2000 
spatial frequency (Hnas/mm) 


Figure 2 Optical transfer function without and with the 
shading band. 


4. Conclusions 

We have shown how a simple shading band modifies the optical transfer function of a data storage 
drive. This could improve contrast in the data signal. With the shading band, both laser and medium 
noise levels are reduced. With the particular shading band we choose, the CNR is reduced at low 
frequencies, unchanged for mid-range frequencies, and improved for high frequencies. Due to the property 
of reducing the noise floor, further improvements in CNR might be realized by a more clever design of the 
shading band. 

5. References 

1. T. Milster and C. Curtis, Applied Optics, 31, pp. 6272-6279 (1992). 



Figure 3. A) electronic noise floor, B) laser + 
electronic noise floor, and C) medium + laser + 
electronic noise floors for our system with and 
without the shading band. 



Figure 4 . carrier, total noise, and narrow-band 
carrier-to-noise (CNR) measurements for our 
system. Note that, with the shading band, the CNR 
is reduced at low frequencies, unchanged for mid 
frequencies, and improved at high frequencies. 
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AB STRA CT 

We show how FWHM, FWl/e 2 , Strehl ratio, and encircled energy figures of merit vary with 
different types of aberration and measurement methods. We examine in detail the array sampling method 
and the slit-scan method. Our irradiance in the exit pupil of the optical system is a simple gaussian. We 
found that in general the slit-scan method and the array method do not yield the same result. The width 
w measurements for the central lobe of the diffraction pattern are very insensitive to aberration. 

2. INTRODUCTION 

As has been known since opticists started making lenses, aberrations redistribute energy from the 
central lobe of the Airy pattern into the sidelobes. However, the study of the effects of aberrations in 
^ laser beam systems with nonuniform irradiance profiles is a relatively new topic [1][2][3]. For 

example, a full-width at half maximum (FWHM) of the central lobe in a focused beam that is near the 
value obtained with an unaberrated system is often used as justification that the beam is "diffraction 
limited.” In this paper, we show that caution must be used in drawing this conclusion because of the 
effects third-order aberrations have on the FWHM. We also describe how several figures of merit vary 
with aberrations and measurement technique. 

Two methods used to measure laser beam quality are sampling 
an array of data points and measuring with a scanning slit [3]. We 
can obtain a two-dimensional array of sample points by inserting a 
CCD camera or a scanning pinhole directly into the beam. Line 
profile and edge response data are obtained by scanning a slit across 
s? the beam. Instruments used to obtain array and scanning slit data are 

typically much simpler and less expensive than interferometers, which 
is a third alternative to measure beam quality. In this paper, we 
concentrate on the properties of array and slit-scan measurements of 
a focused laser beam. 

Common figures of merit for laser beam quality include width 
measurements, Strehl ratio, and encircled energy, as shown in Figure 
1. FWHM and full-width at 1/e 2 (FWl/e 2 ) are self explanatory. The 
Strehl ratio is defined as the ratio of the peak irradiance of the 
aberrated beam, I*, to the peak irradiance, I, of a system with no 
aberrations. Encircled energy is defined as the ratio of the power, P’, 
w contained in a small circular region around the peak of the aberrated 

beam to the total power, P, in the beam. Our circular region is equal 
to the diameter of the first ring of the Airy pattern in an aberration- 
free and uniformly illuminated system. This diameter is 1.22 X/NA, 
where X is the laser wavelength and NA is the numerical aperture of 
the focused beam. A related figure of merit is the energy ratio of the power outside the circular area to 
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laser beam 


the total power in the beam, which can be calculated from: energy ratio = 1 - P*/P. 

Our optical system model is shown in 
Figure 2. A simple gaussian laser beam is 
reimaged through an optical system so that the 
waist is located a distance R from the exit pupil. 

The laser beam irradiance in the exit pupil has a 
FWl/e 2 of 2w. The stop has diameter d. We 
assume that d > > X and that the NA is large 
enough so that focus shifts due to the properties 
of the gaussian beam are insignificant. We also 
assume that the NA is small enough so that vector 
diffraction effects are not significant. We model 
the effects of aberrations by adding a phase error, 

AW, to the wavefront in the exit pupil. The 
phase errors take the form of astigmatism (W^, 
spherical (W^o) and coma (W 131 ). 



Figure 2. Optical system model. 


In the following paragraphs we review array and slit-scan methods. Then we discuss effects of 
individual third-order aberrations on various figures of merit. Next we discuss effects of random 
combinations of third-order aberrations. We then summarize and present our conclusions. 

3. REVIEW OF ARRAY AND SLIT-SCAN METHODS 


A :sDtscan 


B :edge scan 


laser beam 




L 


Array methods include scanning a pinhole 
over the measurement plane or using a two- 
dimensional array of detectors, like those found in 
a charge-coupled device (CCD) camera. In either 
case, a two-dimensional data set is acquired from 
which the figures of merit are calculated. The 
data are discrete samples of the irradiance incident 
onto the measurement plane. In order to provide 
sufficient sampling, the spacing between detector 
elements must be small compared to the beam 
size. Typically, CCD cameras have pixel sizes 
on the order of 10 fim on a side. This limits the 
practical beam sizes that can be measured to 
several hundred microns. An auxiliary lens may 
be used to produce a magnified image of a 

smaller beam onto the detector plane, but additional aberrations are often introduced that affect the 
measurement. Also, care must be taken in interpreting metric information from the array, because pixels 
in CCD cameras are often not square nor do they have the same interval in the horizontal and vertical 
directions. In our study we assume that the pixels are square and uniformly spaced. 

Slit-scan methods are used to derive one-dimensional information from the laser beam. As shown 
in Figure 3A, a narrow slit is used to scan the measurement plane in the x direction. The slit integrates 
the irradiance in the y direction, so data do not represent true beam profiles, but rather they represent 
projections along the y axis. The width of the slit determines the resolution of the measurement. Ideally, 
an infinitely narrow slit would be used, but as the slit becomes too narrow the signal-to-noise degrades. 


Figure 3. 
method. 


(a) slit-scan method, (b) edge scan 
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Another way to obtain slit-scan data is from a moving knife edge, as shown in Figure 3B. In this case, 
an opaque surface with a sharp edge is scanned across the measurement plane in the x direction. The data 
from the edge scan are differentiated, and the result is equivalent to an infinitely narrow slit scan, except 
that the signal-to-noise is improved because more signal light is available. 

4. EFFECTS OF INDIVIDUAL THIRD-ORDER ABERRATIONS 



array and slit-scan 
astigmatism. 


methods as a function of 


In this section we examine the effects of 
astigmatism, spherical aberration and coma on 
several figures of merit. Both array data and slit- 
scan methods are implemented. In our study, we 
use a simple scalar diffraction modeling code. 

The optical system is such that 2w/d = 0.89, 
which corresponds to the optimum overfill in 
terms of maximum peak irradiance in the focused 
beam [4]. In a real optical system, defocus 
and tilt can be adjusted to give a higher beam 
quality by canceling some of the effects of 
aberrations. In our modeling we add the 
appropriate amount of defocus and tilt to 
minimize the root-mean-square (rms) wavefront 
error. The amount of phase error, AW, added to 
the ideal wavefront for each aberration is 
described by the peak error at the edge of the exit 
pupil before correction. Width measurements are normalized with respect to 

NA and X. For a specific optical system, the physical width is found by multiplying the normalized value 
byX/NA. 

We consider the FWHM figure of merit 
first. Figure 4 displays the normalized FWHM 
versus astigmatism for the array and slit-scan 
methods. The laser beam exhibits a symmetric 
profile in the measurement plane because defocus 
has been added to minimize the rms wavefront 
error. There is little change in the FWHM 
measurement for less than one wave of 
astigmatism. The array method and the slit-scan 
method produce similar results. Figure 5 displays 
the normalized FWHM versus spherical 
aberration. Sidelobes due to more than two 
waves of spherical aberration affect width 
measurements significantly. Lobe widths plotted 
in Figure 5 include the maximum FWHM (that 
including all of the sidelobes), central-lobe 
FWHM and side-lobe FWHM. As the amount of aberration increases, the FWHM of the central lobe 
actually decreases for the array method. The slit-scan method is slightly more sensitive. The jump in 
the FWHM around 2.4 waves is due to the shape of the sidelobes in the slit-scan method. As shown in 
Figure 6A, the slit scan of a beam having 2.4 waves of spherical aberration has sidelobes that increase 



Figure 5. Normalized FWHM of the array and slit- 
scan methods as a function of spherical aberration. 
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the FWHM. With 2.3 waves of spherical aberration, the sidelobes are below the half-maximum of the 
irradiance peak. In Figure 6B, we show the array profile of the beam with 2.4 waves of spherical 
aberration. Note that the central lobe is well defined, and the sidelobes are well below the half-maximum 
of the peak irradiance. Figure 7 displays the normalized FWHM versus coma. Coma is not a symmetric 
aberration, so we consider profiles and slit-scans in both the x and y directions. Again, FWHM is 
insensitive below one wave of aberration. Sidelobes affect the FWHM measurement at about 2.6 waves 
of coma. An oscillatory property is observed for the FWHM of the central peak in the x direction, which 
is also the direction of the maximum wavefront error in the uncorrected system. 
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Figure 6a. Irradiance profile of the slit-scan Figure 6b. Irradiance profile of the array 
method with 2.4 waves of spherical aberration. method with 2.4 waves of spherical aberration. 



FWl/e 2 figure of merit becomes complicated to 


interpret above one wave of coma. In general, the slit-scan method is more sensitive. Note that for all 
three aberrations the FWl/e 2 figure of merit is not sensitive below 0.5 wave of aberration. 

We now discuss Strehl ratio. Figure 10 displays the Strehl ratio for all three aberrations and both 
measurement methods. For the array method , Strehl ratio is a very sensitive figure of 
merit. For the slit-scan method, Strehl ratio is not as sensitive, but it does provide a monotonically 
decreasing figure of merit with increased aberration. Both methods, Strehl ratio is most sensitive to 
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astigmatism. 



Figure 8. Normalized FW 1/e 2 of the array and 
slit-scan methods as a function of spherical 
aberration. 



Figure 10. Strehl ratios of the array and slit- 
scan methods as a function of waves of aberration. 



Figure 9. Normalized FW 1/e 2 of the array and 
the slit-scan method as a function of coma. 



Figure 11. Energy ratio of the power outside 
of the central lobe to the total power in the beam. 


Figure 1 1 displays the energy ratio for all three aberrations and both measurement methods. Like 
we observed in the Strehl ratio calculations, the energy ratio is a smooth monotonic function versus 
aberration. The array method is more sensitive than the slit-scan method. Coma and astigmatism are 
more easily detected than spherical aberration. When data in Figure 11 is compared to the width data 
in Figures 4, 5, and 7-9, we observe that the insensitivity of the FWHM can be a severe problem. For 
example, up to 40% of the total spot energy is contained outside the central lobe for one wave of 
aberration. Although FWl/e 2 is slightly more sensitive, up to 18% of the total spot energy is contained 
outside of the central lobe for one-half wave of aberration. 


5, E FFE CT S OF CO MB INED ABERRAT ION S 

We study the effects of combined aberrations with a simple extension to our basic model. Instead 
of a single aberration, we included random amounts of astigmatism, spherical, and coma. In addition. 
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the coma rotation angle and the astigmatism 
rotation angle were included as random variables. 

Tilt and defocus were added in order to minimize 
the nns wavefront departure. Standard deviation 
of the wavefront, a , was kept between the limits 
: 0 £ a £ 0. 10. The overfill of the laser beam 
in the exit pupil was slightly different in the x and 
y directions. A total of 2000 trials were 
performed. The figures of merit for both the 
array method and the slit-scan method were 
calculated for each trial. The result of the 
average FWHM versus a is shown in Figure 12, 
which displays similar characteristics to Figures 5 
and 7 for low values of aberration. The 
difference between the focused beam width in the 
x and y directions is due to the difference in 
overfill in the exit pupil. The slit-scan method 

produces smaller beam widths than the array method for a < 0.07 and larger beam widths for a > 0.08. 
The relation between a and the peak AW depends on the type of aberration present. For example, one 
wave of astigmatism has a — 0.200, one wave of spherical has a = 0.074, and one wave of coma has 
a = 0.116. 



Wavefront Standard Deviation (Waves) 
Figure 12. Normalized FWHM of the array and 
slit-scan methods as a function of wavefront standard 
deviation. 


6 SUMMARY AND CONCLUSIONS 

We have shown how FWHM, FWl/e 2 , Strehl ratio, and encircled energy figures of merit vary 
with different types of aberration and measurement methods. We found that in general the slit-scan 
method and the array method do not yield the same result because the slit-scan method measures an 
integrated line profile of the beam while the array method measures a profile. The FWHM and FWl/e 2 
values for the central lobe of the diffraction pattern are very insensitive to aberration. Therefore, one 
should use caution when claiming that an optical system is "diffraction limited" based solely on the these 
criteria. In the case of spherical aberration the central lobe width actually decreases with increased 
aberration. The slit-scan method width measurements are typically more sensitive to aberration than the 
array method. The most sensitive figure of merit is the array method Strehl ratio. The array method 
energy ratio is a useful figure of merit because is it describes the ratio of the power outside the central 
lobe to the total power in the beam. The numeric values for the energy ratio from a slit-scan are always 
lower than the corresponding array method. Both the Strehl ratio and the energy ratio are smooth, 
monotonic functions versus aberration. A computer experiment in which random combinations of 
aberration are added to the exit pupil indicates that the average FWHM is different for slit-scan and array 
methods. Our results for individual aberrations are based on one condition of overfill in the exit pupil. 
Other overfill ratios (2w/d) could yield different results. 
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Abstract 

The design of a near- field scanning optical microscope (NSOM) is described. Applications to 
the study of magneto-optic media and domain detection are discussed. 
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The ability to achieve subwavelength optical resolution for microscopy through the use 
of a mechanically scanned near-field optical probe has recently been reported. 1 This technique 
provides a means of exceeding the far-field diffraction limit associated with the use of an 
objective lens in conventional optical microscopy, and is known as near-field scanning optical 
microscopy (NSOM). The near- field probe may be realized by forming a subwavelength-sized 
aperture in an opaque screen. Such a probe may be implemented in practice by tapering an 
optical fiber down to the desired diameter and coating its exterior (except for the very tip) with 
aluminum. The spatial resolution achieved is dependent on the aperture diameter and the tip-to- 
sample separation, and the resolution is approximately equal to the aperture diameter for 
separations less than one aperture radius. 

NSOM is capable of providing imagery with polarization contrast and is, therefore, a 
potentially useful tool for the study and analysis of magneto-optic materials. NSOM also 
provides a means of performing very high-density MO data storage, and domains as small as 60 
nm have been successfully written and detected. 2 High-resolution lithography on a submicron 
scale has also been demonstrated. 3 

Our NSOM instrument has been designed for the examination of MO thin film samples 
and is illustrated in Figure 1. The instrument consists of three optical subsystems: a 
conventional visual polarization microscope, a servo subsystem for maintaining the fiber probe 
position relative to the sample, and a characterization subsystem for obtaining the MO signal. 
The sample is mounted on a piezo-driven x-y flexure stage, and characterization is performed 
in transmission. The probe tip may be coarsely prealigned to sample features, such as domain 
structures, with the visual subsystem. The visual polarization microscope operates in reflection 
and consists of the lamp, polarizer, beamsplitter, objective, prism, analyzer, and eyepiece. The 
servo subsystem regulates the separation (typically less than 20 nm) of the tip and sample with 
a hybrid shear atomic force/optical technique. A small dither is applied to the fiber tip (parallel 
to the sample surface), and light from the HeNe laser is scattered from the tip and imaged to the 
pinhole. The dither produces a modulated signal on the PMT. Atomic forces between the tip 
and sample perturb the amplitude and phase of the dither, which can be detected and used as a 
feedback signal to regulate tip position. The characterization subsystem operates at blue 
wavelengths and generates an MO signal in transmission. Light from the Ar laser is launched 
into the fiber probe with the aid of waveplates to produce any arbitrary state of incident 
polarization. The beamsplitter and prism are removed from the visual subsystem, and 
polarization-sensitive optics detect the signal light as the sample is scanned to create an image. 



w Figure 1. NSOM system layout for detection of MO domain structures. 

A typical differential detection scheme is shown consisting of waveplates, a polarizing 
beamsplitter, two PMTs, and a difference amplifier, but other schemes are possible as well. 

_ The samples used in this study are Co/Pt multilayers prepared by sputtering. The design 

w is 13 bilayers of 3.S A Co and 8 A Pt, yielding a total MO film thickness of 149.3 A, on a 

Pyrex substrate. No cover layer is used because the probe tip must be in proximity to the MO 
layer surface to obtain the resolution afforded by the near-field interaction. Domains were 
prewritten on the samples by conventional far-field means using an Olympus polarization 
microscope fitted with a focused laser diode source and an electromagnet. 

w 

NSOM may be applied to a variety of studies of MO materials and data storage 
problems. The high spatial resolution of NSOM (down to 12 nm) can be utilized to study 
domain wall jaggedness. Depolarization effects due to disk grooves and surface roughness can 
be examined on a subwavelength scale. The interaction of the near-field probe and MO domain 
^ structures can be characterized with various detection schemes. Domain-sensitive methods such 

as differential detection with linearly polarized light, and edge-sensitive methods using split 
detectors and circularly polarized light are two techniques of particular interest. 
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: w 


We present a thorough study of differential wax-wane focus servo technique including 
effects of aberration and cancellation of crosstalk. 
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A Complete Study of Differential Wax-wane Focus Servo Technique 

Mark S. Wang and Tom D. Milster 
Optical Sciences Center, University of Arizona 
Tucson, Arizona 85721, (602) 621-8280 

Summary 

We concentrate on differential wax-wane focus servo technique 1 that is insensitive to 
many of the errors found in other schemes. Servo optics are illustrated in Figure 1. The beam 
reflected from optical disk is focused by the detector lens onto two detectors, detector 1 is 
slightly inside focus, and detector 2 is slightly beyond focus. The quad detector is offset from 
the center of the beam to give the FES algebra, where a is an electronic gain factor that can be 
any number larger than one. 

We used scalar diffraction 
model to study its performance. 

We studied beam propagation in 
the optical system, the focus error 
signal, the detector alignment 
tolerance, the tracking error signal 
and crosstalk. The effect of 
aberration on the above 
parameters is modeled in detailed. 

The differential wax-wane 
technique has several advantages 
over a single wax-wane focus 
servo technique. The gain is two 
times higher, the lock-on-range is 
better defined, and the linearity is 
ten times better in terms of RMS 
deviation. It is also insensitive to track rotation and disk tilt. 

One important performance parameter is sensitivity to pattern noise , which is a false 
focus error caused by small changes in the spatial distribution of optical power on the detectors. 
The most common type of pattern noise is due to diffraction from tracks. It is also referred as 
tracking crosstalk. 

The effect of aberration varies depending on where the aberration is introduced and what 
kind of aberration it is. The focus offset in differential channel is the same as in a single 
channel. Crosstalk, change of gain, and nonlinearity are improved in the differential channel. 
W 222 ( 1 ) and W, 3 ,( 1 ) generate no focus offset and don’t change the gain or linearity, but they 
change the amplitude and phase of track crosstalk significantly. W 13 ,( || ) has minimal effect on 
any parameter if it is disk aberration, but changes focus offset, gain, and crosstalk amplitude if 
it is source or detection-optics aberration. The most significant effect of W^f [| ) and is 
the focus offset. 

The case of combined aberrations (in which we measured the aberrations on our optical 
system) generates a large focus offset that is a strong function of the electronic gain factor a. 
There is about 0.7 ^m of track crosstalk in a single channel, and the amplitude and phase are 
different between the two individual channels. Thus, crosstalk is not canceled in differential 



technique, a is an electronic gain factor. 
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channel. The crosstalk amplitude 
of individual channels differs by a 
factor of about 2, and the relative 
phase is also different. Total 
cancellation is not achieved in the 
differential FES. However, we 
found that, by rotating the 
detector, the phase of crosstalk 
can be adjusted so that the 
individual crosstalks are 90° out of 
phase with the tracking signal. 
T , crosstalk phase of detector 1 
is then the same as detector 2. 
We can also change the electronic 
gain to compensate for the 
amplitude differences. By 
reducing the electronic gain of 
detector 1 by a factor of 2, the 
crosstalk amplitude become the 
same. The crosstalk in the 
differential channel is canceled 
very well. Theoretically, the 
crosstalk can be reduced to zero 
by fine tuning the electronic gain 
factor and rotating the quad 
detectors. The residual crosstalk 
is less than 0.1 jim. Figure 3 
shows the measured crosstalk. 
Note that, by changing the 
electronic gain, the differential 
FES is e = e,/2 — e 2 . 
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Figure 2 Crosstalk amplitude as a function of the 
electronic gain factor a when various aberrations 
are present in the optical system. 
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Measured FES crosstalk and TES as a function 
of track shift, t = - e 2 . 


Conclusion: We presented a 
complete study of the differential 
wax-wane focus servo technique 
including both scalar diffraction 
modeling and experimental 
measurements. The crosstalk was 

reduced from 0.7 /xm in a single channel to less than 0.1 fim in the optimized differential 
channel. This gives the differential wax-wane technique a clear advantage over several 
techniques we compared. 
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Introduction: 


Just prior to its being EFM modulated (i.e., converted to eight-to-fourteen channel data by 
the EFM encoder) and written to a Compact Disc (CD), information that passes through 
the CIRC Block Encoder is grouped into 33-byte blocks referred to as EFM frames. 24 of 
the bytes that make up a given EFM frame are user data that was input into the CIRC 
encoder at various (different) times, 4 of the bytes of this same EFM frame were created 
by the C2 ECC encoder (each at a different time), and another 4 were created by the C 1 
ECC encoder (again, each at a different time). The one remaining byte of the given EFM 
frame, which is known as the EFM frame C&D (for Control & Display) byte, carries 
information that identifies which portion of the current disc program track the given EFM 
frame belongs to and also specifies the location of the given EFM frame on the disc (in 
terms of a time stamp that has a resolution of l/75th second, or 98 EFM frames). (Note: 
since the program track and time information is stored as a 98-byte word, a logical group 
consisting of 98 consecutive EFM frames must be read, and their respective C&D bytes 
must be catenated and decoded, before the program track identification and time position 
information that pertains to the entire block of 98 EFM frames can be obtained.) The 
C&D byte is put at the start (Oth byte) of an EFM frame in real time; its placement 
completes the construction of the EFM frame - it is assigned just before the EFM frame 
enters the EFM encoder. In Fig. 1 ., the 33 x 1 column vector EFM <,:> represents EFM 
frame No. "i"; the byte labeled as F 0 j in this column vector is the C&D byte of EFM 
frame No. "i". 

In the 134 (row) x 33 (column) array that accompanies this text, each row represents an 
EFM frame and each EFM frame's identification number is given by the value of the first 
element of each row, i.e., these values are found in the column identified as F 0 at the top 
of the 134 x 33 array. The numerical value of these (column F 0 ) array elements 
corresponds to the value of the second subscript ("i") of the element F 0 j in the column 
vector EFM <i> of Fig. 1. EFM frames are written sequentially to the CD disc when it is 
recorded. Thus, byte F 0 of EFM frame No. 0 is recorded first, followed by byte Dj of 
EFM frame No. 0, ..., and the 33rd byte written to the disc is byte CI 3 of EFM frame No. 
0. As EFM frame No. 1 directly follows EFM frame No. 0, the 34th byte written to the 
disc is byte F 0 of EFM frame No. 1 , etc. The 134 x 33 array is therefore a map of how 
data is written to the disc; the first row (i.e., the top row) of the array is written first, left 
to right, followed by the next row, again left to right, etc. 

In the following discussion we will refer to four distinct blocks of data: (i) 24-byte User 
Input Data blocks, (ii) 28-byte C2 words, (iii) 32-byte Cl words and (iv) 33-byte EFM 
frames. In order to track the location on the recorded CD disc of the various data bytes, 
we shall "tag" (or identify) each data byte that is recorded on the disc with a number, or 
index, that corresponds either to the instant in time at which the particular byte entered 
the CIRC encoder (in the case of User Input Data), or to the instant in time at which the 
data byte was created by the CIRC encoder (as is the case with C2 parity bytes, C 1 parity 
bytes and EFM frame C&D bytes). Once a byte is "tagged", the number it is tagged with 
will be kept, i.e., that number will not change, even when the particular byte becomes a 
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member of other CIRC information blocks (which will occur as the byte flows through 
the CIRC encoder). In order to assign numerical values to the byte "tags", we shall think 
of information flowing through the CIRC encoder in expanding-size blocks, namely, the 
User Input Data blocks, C2 and Cl words, and EFM frames mentioned above. 

The User Input Data Blocks: 

24 bytes of sequential user data is input to the CIRC encoder as a single 24-byte User 
Input Data block. A User Input Data block that enters the CIRC encoder at instant "k" is 
represented by the column vector D <k> of Fig. 1. The sequence in which user data bytes 
filled the User Input Data block is indicated by the first subscript of each data byte (i.e., 
element) of D <k> ; the instant at which the block entered the CIRC encoder is denoted by 
the second subscript of the various bytes of D <k> . Thus, Dj k denotes the first byte of the 
24-byte User Input Data block that entered the CIRC encoder at instant "k" and D 2 4 k 
denotes the 24th (and last) byte of the User Input Data block that entered the CIRC 
encoder at instant "k". Sometimes we shall refer to blocks that enter the CIRC encoder at 
instant "k" as block No. "k", and vice versa. 

The C2 Words: 


Directly after they enter the CIRC encoder, the bytes of a 24-byte User Input Data block 

are reordered (or swapped) such that the sequence (D| k , D 2k , D 3k , , D 23 k , D 24 k } 

becomes the sequence {D l k , D 2 k , D 9 k , D 10 k , D 17 k , b l8k ,’D 3)k ,’D 4k , Dj 1>k , D 12 k , 

D 19,k> D 20,k> D 5,k’ D 6,k> D 13,k> D 14,k> D 21,k> D 22,k> D 7,k> D 8,b D 15,k> D 16,k> D 23,k> D 24,k)- 


Then, after delaying the first 12 bytes of the reordered User Input Data block by 48 bytes 
(which amounts to interleaving the first half of the reordered User Input Data block by 
two reordered User Input Data blocks) the resultant (reordered, two-block interleaved) 
24-byte block is encoded by the C2 ECC coder. This encoding process involves the 
computation of four C2 parity bytes and the subsequent addition of them to the reordered 
and doubly-interleaved 24-byte block between bytes D 2 q k _ 2 and D5 k . This results in a 
28-byte C2 word such as the one represented by the column vector C2 < P > in Fig. 1 . Four 
new bytes, namely the C2 parity bytes C2p 0 , C2pj , C2p 2 , and C2p 3 are created by the 
CIRC encoder when the C2 word is formed; we therefore "tag" these four parity bytes 
with a numerical value that corresponds to the instant at which these parity bytes were 
created. In addition, we "tag" the C2 word itself with this same numerical value. This is 
illustrated in Fig. 1 , where it is shown that the column vector C2 < P > , which represents C2 
word No. "p", contains the four parity bytes C2p 0 p, C2 p ]p , C2p 2 p and C2p 3 p (recall our 
convention is that the second byte subscript contains the numerical value that corresponds 
to the instant at which the C2 word is created, or equivalently, to the C2 word's No. 
"value"). This same C2 word (i.e., C2 word No. "p") also contains bytes from User Input 
Data block No. "p-2" as its first 12 information bytes (these 12 bytes experienced a delay 
of 2 User Input Data blocks prior to entering the C2 ECC encoder), as well as bytes from 
User Input Data block No. "p" as its last 12 information bytes (these 12 bytes experienced 
no delay between the time they entered the CIRC encoder and the time at which they 
entered the C2 encoder - in effect these 12 bytes were input directly into the C2 encoder). 
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The Cl Words: 


Next, the 28-byte block output by th(TC2 coder is interleaved again, this time on a byte- 
by-byte basis. Here, the 2nd byte of a given C2 word is delayed by four C2 words (112 
bytes) relative to the first byte, the 3rd byte of the given word is delayed by eight C2 
words (224 bytes) relative to the first byte, ..., and the 28th byte of the given word is 
delayed by one-hundred-and-eight C2 words (3,024 bytes) relative to the first byte. A 28- 
byte block of data that has been variably-interleaved in this manner is input to the C 1 
ECC encoder, which computes four Cl parity bytes and adds them to the end of the 
variably-interleaved 28-byte block that was input to it. This is illustrated in Fig. 1, where 
it is shown that column vector Cl <n> , which represents Cl word No. "n", has Cl parity 
bytes Clpo n , Clp ]n , Clp 2>n and Clp 3p as its last four bytes. From Fig. 1 we also see 
that the 28 information bytes of Cl word No. "n" (i.e., its first 28 bytes) consist of user 
data bytes taken from User Input Data blocks that entered the CIRC encoder at various 
times ranging from two Cl word-times prior to the time of the current word (cfi, the first 
byte of Cl <n> , namely D, n _ 2 , which entered the CIRC encoder as a member of User 
Input Data block No. "n-2") to one-hundred-and-eight Cl word-times prior to the current 
word (cfi, the 28th byte of Cl <n> , namely D 24n _ 108 , which entered the CIRC encoder as a 
member of User Input Data block No. "n-108"). It can also be seen from Fig. 1 that the 
information portion of a Cl word contains parity bytes that belong to four different C2 
words (namely, bytes C2po n . 4 g through C2p 3n _6o of column vector Cl <n> ). 

The EFM frames: 

A final interleaving follows the Cl encoding. The odd bytes, namely the 1st, 3rd, 5th, ..., 
31st bytes, of the 32-byte block that is output by the Cl coder are each delayed by one Cl 
word (i.e., by 32-bytes) relative to the even bytes of the same Cl word. And lastly, a 
single EFM C&D byte is appended to the beginning of each of the resultant, singly- 
interleaved 32-byte blocks to form the 3 3 -byte EFM frames. The column vector EFM <i:> 
in Fig. 1 is a representation of EFM frame No. "i". Note that this EFM frame's C&D 
byte, namely byte F 0 i? and the two Cl parity bytes C 1 pj j and Clp 3 i , are the only bytes 
in EFM frame No. "i" that were created (by the CIRC encoder) at instant "i". The EFM 
frame C&D byte Fq , is created by the CIRC encoder at the instant the EFM frame No. "i" 
is formed, and the two Cl parity bytes mentioned directly above, were created by the C 1 
encoder at this same time, which is why they are labeled as parity bytes that belong to Cl 
word No. "i". 

Locating Specific Information Bytes and CIRC Data Structures on the Recorded 
Disc: 

In order to keep track of the various processing steps outlined above, we can use an index 
that increments once each time a new EFM frame is created. Concurrent with that event 
(i.e., the creation of an EFM frame), a new 24-byte User Input Data block is input to the 
CIRC encoding system, a new C2 word is encoded by the C2 coder, and a new Cl word 
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is encoded by the Cl coder. Thus, we tag all 24 bytes of a single User Input Data block 
with a common index, the four C2 parity bytes that correspond to a given C2 word share 
a common index (since they were all created at a particular time), and the four C 1 parities 
of a given Cl word are given a common index (again, since they were simultaneously 
created). Specifically, the 24 bytes of a given User Input Data block are all given the 
index value assigned to the EFM frame C&D byte that was current when the first byte of 
the 24-byte User Input Data block entered the CIRC encoding system. Similarly, the 
group of four C2 parity bytes, and the group of four Cl parities that were each created via 
single encoding operations, are each given the index value of the EFM frame C&D byte 
that was current when those groups of parities were computed and assigned by the 
respective ECC coders. 

It is these index values that are listed in the 134 x 33 array, which is our map of the 
recorded CD disc. From that array we see that the C&D byte of EFM frame No. 0 (byte 
Fo in the first, i.e., top, row of the array) is assigned an index = 0, that user data byte #1 of 
EFM frame No. 0 (byte D] in row 1) came from the User Input Data block having index = 
-3 (i.e., this byte came from the User Input Data block that was input while the EFM 
frame with index = -3 was current), that user data byte #2 (byte D2 in row 1) came from 
the User Input Data block having index = -6 (i.e., the User Input Data block that was 
input while the EFM frame with index = -6 was current), that user data byte #9 (byte D 9 
in row 1) came from the User Input Data block having index = -11, ..., that C2 parity byte 
number 0 (byte C2 q in row 1 ) was created by the C2 coder when the EFM frame with 
index = -49 was current, ..., and that Cl parity byte number 3 (byte CI3 in row 1) was 
created by the C 1 coder when the EFM frame with index = 0 was current. 

The data in a given EFM frame is sequentially recorded on the disc and sequentially 
indexed EFM frames are recorded in sequence. Thus, the 134 x 33 array is essentially a 
map of how data is recorded (laid out) on the disc. Byte Fo of EFM frame No. 0 is 
recorded first, followed by byte Dj of EFM frame No. 0, etc. Byte Cl 3 of EFM frame 
No. 0 is the last byte of that frame to be recorded and it is directly followed by byte F 0 of 
EFM frame No. 1 . Because of this, the array can be used to determine the location of all 
the data from a specific User Input Data block, or from a specific C2 or Cl word. For 
example, the 24 data bytes located in columns 2 through 1 3 and in columns 1 8 through 29 
that have the same index value all came from the same User Input Data block (e.g., byte 
D) in EFM frame No. 0, byte D 2 in EFM frame No. 3 and byte D 9 in EFM frame No. 8 
all are members of the User Input Data block having index = -3). With the help of the 
column vectors given in Fig. 1, the locations on the disc of all the bytes of specific C2 
and C 1 ECC words can also be determined. 

Using the column vector representations for an arbitrary EFM frame, an arbitrary C2 
word and an arbitrary Cl word that are given in Fig. 1, we see that C2 word No. "p" (i.e., 
the C2 word that was created when EFM frame No. "p" was current) has data bytes that 
were input as part of User Input Data block No. "p-2" as its first 12 bytes and data bytes 
that were input as part of User Input Data block No. "p" as its last 12 bytes. Thus, if one 
wanted to read the disc and retrieve, for example, all the data (24 bytes) that were input as 



User Input Data block No. 0, one would have to read a length of disc track which 
contains the entire 56 bytes that make up C2 words No. 0 and No. 2 (since the last 12 
bytes of C2 word No. 0 and the first 12 bytes of C2 word No. 2 contain the 24 bytes that 
were input as User Input Data block No. 0). To acquire C2 word No. 0 one would have 
to read all those bytes in columns 2-13 (bytes D] through D20) of the 134 x 33 array that 
have index = -2 and all those bytes in columns 14-29 (bytes C2 0 through D24) that have 
index = 0. Similarly, to acquire C2 word No. 2, one must read all those bytes in columns 
2-13 that have index = 0 and all those bytes in columns 14-29 that have index = 2. 

From the 134x33 array we see that this means that EFM frames No. 1 through No. 1 1 0 
must be read. However, note that byte D24 of EFM frame No. 1 10 (which is the last byte 
of C2 word No. 2) is also part of C 1 word No. 1 1 0, which has all of its odd-numbered 
bytes in EFM frame No. 1 1 1 (the column vector representation of an arbitrary Cl word 
given in Fig. 1 shows that the 2nd, 4th, 6th, ..., 32nd bytes of Cl word No. "n", i.e., bytes 
D 2,n-6’ D l0,n-14> D l8,n-22> cl P2,n ^ of EFM frame No. "n", while the 1st, 3rd, 
5th, ..., 31st bytes of Cl word "n", i.e., bytes Dj n . 2 , D 9 n _ 10 , D| 7n . 18 , ..., Clp 3 n are part 
of EFM frame No. "n+1 "). This means that EFM frame No. 1 1 1 must be read too, so that 
Cl word No. 1 10 can be acquired and decoded (because Cl decoding must be 
accomplished before C2 decoding). In total, then, 1 1 1 EFM frames must be read from the 
disc in order to retrieve the 24 data bytes that were input as a single User Input Data 
block... 

To augment the understanding of how to use the 134 x 33 array to determine the locations 
of CIRC -encoding data structures on a recorded CD disc, on that array we have drawn 
circles around the locations of all 28 bytes of C2 word No. 0, and drawn squares around 
the locations of all 32 bytes of Cl word No. 110. Triangles have also been drawn around 
the locations of all 24 bytes of User Input Data block No.5. (These CIRC information 
blocks were arbitrarily chosen.) 

CD Audio Interpolation: 

The above analysis sheds light on the purpose for performing the byte reordering 
(swapping) that is done on the Input User Data block prior to C2 encoding. The CD 
system is designed, at its lowest level, as a serial (tape-like) digital audio recorder - audio 
program material is sampled at a rate of 44.1 KHz and those samples are digitized as 32- 
bit words (each 4-byte digital audio sample consists of a 16-bit right-channel audel and a 
16-bit left-channel audel). Thus, a given 24-byte User Input Data block contains exactly 
6 audio samples; the first four bytes (bytes Dj ; - D 4 j, in the case of Input User frame No. 
"j") constitute the first audio sample, the next four bytes (D 3 j - D g j) make up the second 
audio sample, etc. The column vector representation of a C2 word in Fig. 1 shows that 
the first 12 bytes of C2 word No. "p" (which consist of bytes Dj p _2 to D4 p _2, D 9 p _2 to 
Dj 2 p-? and D17 p_2 to D20 p -2 > n scrambled form) would carry the 1 st, 3rd and 5th audio 
samples from User Input Data block No. "p-2", while the last 12 bytes of that same C2 
word (which consist of bytes D5 p to D g p , Dj 3 p to Di6,p and ^21, p t0 ^24, p * n scrambled 
form) would carry the 2nd, 4th and 6th audio samples from User Input Data block No. 

"p". This means that the audio samples carried in a single undecodable C2 word could be 


180 



interpolated from their "neighbor" audio samples, which are carried by two other C2 
words, one of which is recorded at a distance of two C2 words earlier than the non- 
decodable word and the other at a distance of two C2 words later than the non-decodable 
word. For example, if C2 word No. 10 could not be decoded, one could interpolate 
values for audio samples 1, 3, and 5 that belong to User Input Data block No. 8 (and 
which are carried as the first 1 2 information bytes of C2 word No. 1 0) using the audio 
samples 2 , 4, and 6 of User Input Data block No. 8 (which are carried as the last 12 
information bytes of C2 word No. 8 ). Similarly, values for the audio samples 2, 4, and 6 
from User Input Data block No. 10 (which are carried as the last 12 information bytes of 
C2 word No. 10) could be interpolated using the audio samples 1, 3, and 5 from User 
Input Data block No. 10 (which are carried as the first 12 information bytes of C2 word 
No. 12). 

Note that the byte reordering that is done prior to C2 encoding also causes the 
interleaving of the left and right stereo audels that constitute the 1st, 3rd and 5th digital 
audio samples from a given User Input Data block. For example, from the Fig. 1 column 
vector representation of a C2 word, note that the first two bytes (i.e., the left channel 
audel) of digital audio sample 1 , viz., D] and D 2 , are followed by the first two bytes 
(again, the left channel audel) of digital audio sample 3, viz., D 9 and D 10 , which are 
followed by the first two bytes of digital audio sample 5, viz., D 17 and D 18 ; these are 
followed, in sequence, by the last two bytes (i.e., the right channel audel) of digital audio 
samples 1 , 3 and 5 from the same User Input Data block. Digital audio samples 2, 4, and 
6 from a given User Input Data block are treated similarly. This feature allows 
independent (separate) interpolation of values for either left channel or right channel 
audels using "neighbor" left channel or right channel information that is carried in another 
C2 word. 

Handling of Long Burst Errors: 

From the 134 x 33 array of recorded byte index values, we see that a burst error having 
minimum length of 1 00 bytes would be required cause two errors to occur in a single C2 
word. For example, an error burst extending from byte D ( 3 (the second byte of EFM 
frame No. 3) to byte D 2 6 (the third byte of EFM frame No. 6) would contaminate two 
bytes of C2 word No. 2. We can also see that an error burst that has a minimum length of 
1 66 bytes would be required to cause two errors in both of the two C2 words that carry 
"neighbor" audio samples. For example, a burst extending from byte Dj 3 to byte D 2 g 
would be required to cause two errors in both C2 word No. 2 and C2 word No. 4. 

Clearly, very long burst errors can be corrected by the highly-interleaved C2 code, and 
extremely long burst errors are required to defeat the concealment of audio information 
errors via interpolation. 

The C2 code's immunity to long burst errors, however, can readily be compromised if 
there is a high number of short random errors. For example, suppose that a 20-byte long 
burst error occurs; it will cause (only) one byte to be erroneous in each of twenty different 
C2 words. Then, the probability that a short random error will contaminate a second byte 
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of any one of these same C2 words is roughly 20 x 27 x Pg = 540 Pg, where Pg is the 
byte error rate that describes the occurrence of short random errors. Similarly, if a long 
burst error that contaminates two bytes of a given C2 word has occurred, then the 
probability that a short random error will contaminate a third byte of that same C2 word 
is roughly 26 x Pg. Clearly, if the short random error rate at the input to the C2 decoder 
is as high as, say Pg = 1 0" 5 to 1 O' 3 erroneous bytes/byte (a range of rate values that may 
well be representative of the short random error rate of an aged, well-used CD disc), then 
there would be an unacceptably high probability of occurrence of uncorrectable C2 words 
(i.e., C2 words having more than two or three flagged errors), and/or misdecodable C2 
words (i.e., C2 words having more than two unflagged errors). It is the job of the C 1 
code, which is decoded prior to C2 decoding, to reduce the short random error rate at the 
input of the C2 decoder to an acceptable level. 

Correction of Short (Random) Errors: 

The purpose of the single-frame interleaving that is done just prior to Cl encoding is to 
enhance the short random error correcting capability of the Cl code. Specifically, since 
the conditional probability (conditioned on the occurrence of an uncorrectable error 
pattern) of Cl misdecoding, i.e., having the Cl decoder unknowingly output incorrect 
data, is about 1% when a Cl decoder is configured to correct two errors, versus about 
0.0002% when a Cl decoder is configured to correct only one error, it is desirable to use 
single-error correction at the C 1 decoding level. (Since the Hamming distance of both the 
Cl and C2 codes is five, when a codeword of either of these codes is misdecoded the data 
block output by the decoder will have a minimum of five incorrect bytes; thus an 
uncorrectable error pattern that contains, say, only three erroneous bytes will be enlarged 
to a five, or more, byte error if erroneous decoding occurs.) The one-word interleaving 
done at the Cl level allows a single-error-correcting Cl decoder to correct two-byte long 
errors (since such errors will contaminate only one byte in each of two successive C 1 
words). This is important since it has been observed that the average length of the short 
random errors that occur in injection molded CD Digital Audio Discs and CD-ROM discs 
is about 1 1 bits in length. 

The Cl block error rate (BLER) is equal to the fraction of Cl words (at the input to the 
Cl decoder) that contain one or more single or multiple-byte errors. If the occurrence of 
any type of error event is a relatively rare event, e.g., if only 0.1%, or less, of the bytes on 
a disc are erroneous, then the BLER provides a reasonable estimate of the rate of 
occurrence of error events (of any type). However, a long burst error that has length BL 
only increments the BLER count by at most [BL/32], where [x] represents the least 
integer > x. Therefore, the BLER is essentially a measure of the rate of occurance of 
short random errors. 

The Red Book BLER specification is BLER < 0.03. Since there are (at lx playback 
speed) 7,350 Cl blocks/second, this translates to 

BLER < 220 contaminated Cl words/second. 
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In what follows, we shall assume that an error-contaminated EFM frame experiences, on 
average, a single short random error event (i.e., that the short random error rate is 
sufficiently low that the probability of two error events occurring in a single EFM frame 
is very small relative to the probability of occurrence of only one error event). We shall 
also assume that the average length of a short random error is 1.375 bytes (1 1 bits). Such 
an average length error may start at any one of eight locations in an arbitrary byte of a 
depth-one interleaved Cl word. Thus, the 1 .375-byte error event will contaminate only 
one byte in each of the two depth-one interleaved Cl words with probability 5/8 and will 
contaminate two bytes in one, but only one byte in the other, of the pair of interleaved C 1 
words with probability 3/8. 

Random short error rate at the input of the Cl decoder: Since, on average, a random 
short error event always contaminates two Cl words, the short error event rate can be 
approximated by 


BLER short error events 
2 sec 

The coresponding byte error rate is 
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1 events /l erroneous bytes | 

f 1 1 
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1 Cl word 
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p BLER erroneous bytes 
B “ 34,211 byte ' 

Taking BLER < 220 error events/sec, the random short error rate at the input of the Cl 
decoder is found to be 


P B < 6.5 x 1 0~ 4 


erroneous bytes 
byte 


Random short error rate at the output of the Cl decoder: We saw earlier that, on 
average, a short random error event will always cause two Cl words to be contaminated 
and that 3/8 ths of the time one of the two contaminated words will have two bytes in 
error. That is, on average, only one-half of three-eighths of the contaminated Cl words 
will contain more than one erroneous byte. Thus, the byte error rate (due to short random 
errors having an average length of 1 1 bits) at the output of a single-error-correcting Cl 
decoder is given by 



n f 3) BLER \ undecodables ( 2 ^ erroneous bytes / undecodable 

P ’ m lin— J sec X \32 x 7350 J bytes /sec 

When the BLER is at its maximum specified (in the Red Book) value of 220 C 1 
words/second, a single-error correcting Cl decoder will fail to correct about 41 Cl 
words/second. And if a CD disc that is just within the Red Book BLER specification is 
played, the short random error rate at the output of a single-error-correcting Cl decoder is 
estimated to be 


1 


p „ 41x2 

B 1.32x7350 

The short random error rate at the output of a more aggressive Cl decoder, e.g., a double- 
error corrector, would be substantially lower since only those Cl words that have more 
than two bytes in error would remain uncorrected (but, keep in mind that such a decoder 
will have a relatively high probability of misdecoding when an uncorrectable error pattern 
is encountered). Other alternatives for more powerful Cl decoding would be to use the 
Cl decoder to correct a single error together with one erasure per Cl input word 
(assuming, of course, that a means of reliably erasing bytes in the input C l word is 
available), or to use the Cl decoder to correct zero (or one) errors per input Cl word and 
simultaneously detect up to four (two) errors and mark them as erasures for subsequent 
correction by an erasure-correcting C2 decoder. 


erroneous bytes / sec 
bytes / sec 


= 3.5x10" 


erroneous bytes 
byte 


The purpose of the preceding discussion of the short random error rates at both the input 
and output of the Cl decoder was to illustrate how the singly-interleaved Cl code 
contributes to the overall data reliability of the CIRC coding system, and to make the 
point that multiple (> 2) byte errors in C2 words (i.e., at the input to the C2 decoder) 
become highly probable when a long burst error occurs in conjunction with a moderately 
high rate of short random error events (i.e., a moderately high BLER value). 


184 



EFM frame No. V.. 


Cl word No. "n"... 



Figure 1 . The CIRC encoder blocks 
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